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I. INTRODUCTION 


This report describes a numerical procedure permitting the rapid determination 
of the internal performance of a class of scramjet nozzle configurations. The 
geometric complexity of these configurations rules out attempts to employ con- 
ventional nozzle design procedures , wherein properties at the nozzle exit 
plane are specified and wave cancellation techniques are then employed to de- 
sign the wall surfaces. It is not feasible to stipulate exit conditions 
a priori and wave cancellation techniques employing three dimensional character- 
istics are beyond the current state of the art. 

The approach developed is based on the construction of quasi two dimensional 
simple wave networks, wherein lateral expansion effects are incorporated via 
one dimensional approximations, as first suggested by Dr. Antonio Ferri 2 . A 
numerical procedure following this approach has been developed and results 
obtained are highly comparable to those obtained employing a characteristic 
procedure. 

The numerical program developed permits the parametric variation of cowl length, 
turning angles on the cowl and vehicle undersurface and lateral expansion and 
is subject to fixed constraints such as the vehicle length and nozzle exit 
height. The program requires uniform initial conditions at the burner exit 
station and yields the location of all predominant wave zones, accounting for 
lateral expansion effects. In addition, the program yields the detailed pres- 
sure distribution on the cowl and vehicle undersurface and calculates the nozzle 
thrust, lift and pitching moment. 

Due to the differing techniques required for the calculation of frozen flows as 
compared to equilibrium flows, two separate numerical programs have been de- 
veloped. The first program analyzes constant y frozen flow fields and a listing 
of this program is provided in Appendix II. The second program analyzes 
equilibrium hydrogen-air flow fields via equilibrium curve fits and its listing 
is provided in Appendix III. A complete program description is provided in 
Appendix I. The frozen flow program computes a complete nozzle flow field in 
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a fraction of a second while the equilibrium program has a running time of 
several seconds on a CDC 6600, as compared to running times of the order of 

3 

a minute for a quasi two dimensional characteristic approach and twenty 

. 4 

minutes for a three dimensional numerical calculation . 
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II, NUMERICAL PROCEDURES 

Consider a typical nozzle configuration as depicted in Figure (1}, where the 
lateral expansion distribution Z (x) may result from a combination of several 
nozzles merging into a single nozzle. It is assumed in this preliminary 
analysis that the jets after merging are bounded by sidewalls which extend 
downstream of the merged section. The initial flow (at the burner exit) is 
represented as an average uniform flow. The assessment of nonuniformities 
at the entrance station may be obtained applying the numerical procedures de- 
scribed in Reference (3). 

For this configuration, the total amount of expansion from entrance conditions 
is known (based on two dimensional considerations) at the grid points labeled 
V ^ , Cj, A, C, D, Vg, Cgt Eg, E^, F and G and can be readily obtained at 
points Vg, Cg, and E^. The numerical procedure predicts the location of 
these grid points based on three dimensional flow considerations and ef- 
fectively distributes the waves on the cowl and vehicle undersurface to assess 
the pressure distribution. It should be noted that a significant amount of 
the logical procedures in the numerical program are employed to distinguish 
the varying types of wave situations that may be encountered. In the con- 
figurations shown the expansion waves emanating from the cowl (at C^) and the 
vehicle (at V^) are both only partially captured on the undersurface and cowl 
respectively, t , denotes the portion of the cowl expansion wave Av captured 

V v 

on the vehicle undersurface while denotes that portion of the vehicle ex- 
pansion wave Av y captured on the cowl. In other situations, these waves may 
be totally captured or not captured at all, hence all these various situations 
must be distinguishable and treated discretely. 

Before treating the overall logic entailed in a nozzle calculation, the calcula 
tional procedure for a typical grid point will be described: 

A. Frozen Flow Grid Point Calculation - Consider the calcul ational 
procedure required to determine the location and properties of a point 3 (as 
shown in Figure 2) where properties at 1 and 2 are known and 1-3 and 2-3 are 
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characteristic surfaces. The total amount of expansion (Avg) based on two 
dimensional considerations and the flow deflection (e^) is known at point 3. 
Along 1-3 or 2-3 we have 


y 3~ y l,2 

x 3' x l,2 


- g [tan(e 1>2 


± u 1 2 ) + tan(e 3 ± 



(la,b) 


Employing the two dimensional value of expansion Av^ (from initial condition 
i) the Mach. number Mg is obtained via the Prandtl-Meyer relation (where i de- 
notes uniform initial flow properties at the burner exit) 


Av 3 = v 3 -v i 


- (tan ’ 1 ' tan_1 /( ^ (Mr 1 )) 


( 2 ) 


- (tan -1 .tog-l - tan ' 1 41?-1) 


employing an iterative procedure to solve this transcendental equation for Mg.' 
Then, with the Mach angle determined 


y 3 


sin~* 


J 

M 


(3) 


Equations (la) and (lb) yield a tentative location for point 3, and the area 
ratio (A/A*) g is calculated based on two dimensional considerations 




_ 3 


M 3 (1 + M |) 

Y+l 

(I +1 } 


Y+l 

2CT 


(4) 


This ratio is corrected for lateral expansion by multiplying it by the ratio 
Z 3 /Zj r ^here- the‘ lateral expansion variable is expressed by' a suitable poly- 
nomial curve fit 


Z(x) = Ax^ + Bx .+■ C , (5 

where Zg - Z(Xg) and denotes the lateral extent of the nozzle at the initial 
station. 
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(f> = 

* 3 


(f> 




( 6 ) 


Then, the three dimensional corrected Mach number is obtained by replacing the 

two dimensional area ratio in Equation (4) by the three dimensional value 

given by Equation (6), and solving Equation (4) for M* by an iterative 

J3D 

process. Equations (la) and (lb) are resolved using the corrected Mach angle 
W ^3D anCi entl * re P rocec > ure is repeated until two successive values of 
are within a prescribed tolerance. 


t 


A similar procedure is used to determine properties at grid points on boundaries 
with Equation (la) or (lb) replaced with an equation describing the body geom- 
etry. Desired variables (P, T etc.) are then simply obtained by isentropic, 
constant y expansions from initial conditions. 


B. Equilibrium Flow Grid Point Calculation - The geometric location 
of point 3 is obtained employing Equations (la) and (lb) just as for the 
frozen calculation. The known two dimensional value of expansion Av 3 is sub- 
divided into a series of small Av^ increments. The initial value of isen- 
tropic exponent r is obtained from 

r - r (P, h) ' ( 7 ) 

where Equation (7) has been curve fit for equilibrium hydrogen-air as described 
in Reference (5). 


The characteristic compatibility relation 

d In P + dv _ q 
r “ siny cosy 


(8) 


applied across the interval &v. yields the pressure, holding r and p equal to 

J 

their values at the start of the increment. The density is obtained from the 
isentropic pressure-density relation 
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P/p = constant 


(9) 


The velocity is obtained from the Bernoulli relation; 


— + i dV 2 = 0 
P ^ 


( 10 ) 


the enthalpy from the constancy of stagnation enthalpy; 


h + = H » constant 


( 11 ) 


and the Mach number from; 

M = V/a; a = (“) 1 


( 12 ) 


where r has been reevaluated employing Equation (7). This procedure is re- 
peated in small steps Avj until the full wave Av^ has been integrated. 
Having the two dimensional value of point 3 can be tentatively located 
employing Equations (la) and (lb). Then, the two dimensional area ratio 
can be computed from mass flow considerations 

( W 2D - P1VP3Y3 

Since the effective area based on three dimensional considerations is 


(13) 


( A3/A4 ) - (a 3 /aj * z 3 /z 1 

0 1 3D 2D 

the product must be divided by 1 ^/ 1 \ to conserve mass.fl 


ow 


(14) 


(p 3 V,) = (p, V ) * 

43 5 3D J 0 2D j 

Then an iteration procedure is performed to determine the value of three dimen- 
sional expansion (AvO the correct value being that which yields (p~ V 0 ) 

^ 3D ^ ^ 3D 

after application of the integration procedure of Equations (7) thru (12), 

and an update of the location of point 3 using Equations (la) and (lb). 


( 15 ) 
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III. DETERMINATION OF NOZZLE FLOW FIELDS 

A nozzle calculation is performed subject to the following constraints: 

1. The initial profile is uniform. For the frozen flow (con- 
stant y) calculation this requires specification of the pres- 
sure P.. » flow deflection angle e. » Mach number M.. , and spe- 
cific heat ratio y* For the equilibrium calculation one must 
specify P^, e^, M^ , the temperature T.. and the fuel-air 
equivalence ratio 

2. The initial turning at the vehicle undersurface (Av v ) and cowl 
(Av ) occur via sharp corners as depicted in Figure (1). 

3. The wall segments downstream of these sharp corners remain 
straight until the expansion waves emanating from the cowl 
and vehicle undersurface reach the walls (points V 3 and of 
Figure 1). 

4. The nozzle exit height is specified -y_ ). 

v 2 c 2 

5. The recompression on the vehicle undersurface (between and 
V 2 ) is parabolic while the cowl between and C 2 is straight. 

6. The lateral expansion Z(x) is specified via a geometric curve 
fit. 

7. The cowl length and vehicle length are specified. 

The numerical logic employed in the parametric design procedure is to treat 

the cowl length (x -x_ } and the vehicle undersurface expansion Av as para- 

^1 v 

metric variables for fixed values of cowl turning Av c » nozzle exit height, 

lateral expansion and vehicle length. Initially, a short cowl length should be 
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chosen such that the expansion waves from the vehicle expansion fan miss the 
cowl. For this cowl length, the value of the vehicle undersurface expansion 
wave is varied in small increments, the minimum amount of turning being that 
which introduces no recompression in the region V 3 to Vg (i.e., the under- 
surface is flat) to a value for which the recompression produces zero deflec- 
tion at the end of the vehicle. This is illustrated in Figure (3). Then the 
cowl length is increased in specified increments and the entire procedure is 
repeated. 

For a given nozzle configuration, the calculational procedure is as follows: 

1. Grid points A, B, IT and D are calculated following the 
procedure described in the previous section. The two 
dimensional values of expansion from initial conditions 
and flow deflection angles at these points are: 


= 0 , e A = e. 


Av B = Av C' e B ~ e i ' Av C 


y 


AVq“ = AVy, 0-W- = 0.j + AVy 


iv D = Av C + Au V 1 ’ 


= 0^+AVy - A\>£ 


2. The amount of the vehicle expansion wave captured on the cowl (t' c ) 
is determined by an iterative procedure. A value of t‘ c is 
assumed which yields the location of grid points T and W as de- 
picted in Figure (1). The correct value of t c is that for which 
the wave from W intersects the end of the cowl precisely. A 
linear error extrapolation routine is used to speed convergence 
of this iterative process. Having determined t c » the two dimen- 
sional expansion at the end of the cowl (assuming no recompressions 
between and C^) is 
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Av + ?T 
C C 


since the portion of the wave captured is fully reflected. 


3. The amount of the cowl expansion wave captured on the vehicle 
undersurface t v is determined by a similar iterative procedure 
and the two dimensional expansion at the end of the vehicle 
.. is given by 



% + 2t v + Av r 


where the re compress ion Av^ is given by 





Note that for given nozzle parameters, the equation for the para- 
bolic recompression between and 1S determined in the co.urse 
of the numerical calculations. Specification of the cowl length 
and expansion Av yields x„ and y ; the nozzle length and exit 

U Cp Co 

height yields x u and y and the vehicle turn Av yields x and 
*2 Vp v Vo 

y v . Hence, the parabola is determined for each case satisfying 

thi conditions 


y and e specified at x, 
v 3 v 3 J 


specified at x^ 


4 . 


(where 0 = e. + Av ). 

v 3 1 v 

The determination of the exit profile (E Q -E 5 ) is highly dependent 
on the particular configuration being studied. For the case de- 
picted in Figure (1), the cowl and vehicle expansion waves are 
partially captured resulting in a zone of uniform two dimensional 

wave strength between E 9 and E, {Av r - Avr = Av + Av ). Points 

£ j to v c 
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E^ and E^ have this value of expansion and additionally re- 
ceive portions of the reflected waves. Hence, the two di- 
mensional value of expansion experienced at E^ is given by 


Av F - Av + Av + f* t 
L 1 V 1 C 1 c 


where (y P -y r ) 

c 2 

'■‘-Tjpjfr 

t 2 c 2 


while that at E^ is given by 


Av F - Av + Av + g* (t + Av ) 
l 4 V 1 C 1 v r 


where (y F -y v ) 

„ _ , h 4 v 2 
9 " ' 


The corresponding flow deflection angles are: 



6 i 


+ Av w - Av„ - f* 
V 1 C 1 



+ Av y ^ - Av^ + g* (t' v + Av r ) 


Now consider the wave configuration depicted in Figure (4) (where 
the F point and G point are downstream of the exit plane and 
hence not required in the calculational procedure). In this situa 
tion, the two dimensional value of expansion experienced at is 
given by 


= Av Cl + T c + f l" - t c ) 


where f^ 
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while that at is given by 


A \ = A \ + T v + 9 1* (Av c 1 ‘ T v ) 

where y F -y„ 
b 4 ’2 


Si 


y f ~y w 

fc 3 v 2 


The corresponding flow deflection angles are 


e E 1 = e i " Av Cl + T c + f l* < A V T c } 


% = e i + A \ - T v - 9 1* « A V T v>- 


The program developed analyzes both these wave configurations and combinations 
of the two. In addition, it can analyze cases where the full cowl expansion 
is captured on the vehicle undersurface as depicted in Figure (5), and cases 
where no region of uniform two dimensional wave strength on the exit plane 
exists as depicted in Figure (6). For these latter cases, not all points at 
the exit plane are computed by the program due to the logical complexities 
imposed by the wide variety of wave configurations that may occur. 
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IV. DETERMINATION OF THRUST, LIFT AND PITCHING MOMENT 


The thrust, lift and pitching moment are determined by a combination of in- 
tegrating the pressure distribution on the cowl and vehicle undersurface and 
by a momentum flux calculation. Referring to Figure (1), the pressure dis- 
tribution on the vehicle undersurface between Vj and Vg and on the cowl be- 
tween and Cg are readily determined since the two dimensional values of. 
expansion on these surfaces are constant (Av y and Av c respectively). Then, 
the pressure is determined by accounting for lateral expansion by the pro- 
cedure described in Section II (i.e., by employing Equations (1) - (6) for 
frozen flows and Equations (7) - (15) for equilibrium flows with the value 
of two dimensional expansion being constant on each surface). 


Since the lift and pitching moment are quite sensitive to the pressure dis- 
tribution the expansion fans from the cowl and undersurface must be dis- 
tributed properly in the region Vg to Vg on the undersurface and Cg to Cg on 
the cowl. For points on the vehicle undersurface (as depicted in Figure 7) 
the value of two dimensional expansion at a point x* is given by 


* * * 

Av' = Av w + Av^ + 2t 
v r 


where 


and 


* * 

Av =0-8 
r v. 


* -1 
t = tan 


(x -x ) 
v 3 C 1 

(x*-x ) 
C 1 


- 1 . 


(\ -x c ) 

(x v -x ) 
v 2 C 1 


- 1 . 


* tan (t v ) 


(16) 


In the above relation t varies from o to t 1# over the interval x to x , If 

v v 3 m ? 

the full wave is captured on the undersurface as depicted in Figure (5)7 then 
t would vary from o to r, in the interval x„ to x which is affected in 


Equation (16) by replacing x y by x y , 

2 A 


* 

Then t =t. 


between x.. and x 


Similarly, for points on the cowl (Figure 8) the value of two dimensional ex- 
pansion at a point x is given by 



TR 188 


Page 13 


Av = Av c + 2t 


where 


'v , -1 

t = tan 


(x c A ) 

c 3 V 1 


(x-X ) 

1 - V 1 


- 1 . 


{X C A > 
c 3 V 1 


I V r vr r 


- 1 . 


* tan(x c ) 


(17) 


In this relation t varies from o to t in the interval x to x . Knowing 

— C o C p 

the two dimensional values of expansion, the pressure distribution accounting 
for lateral effects is obtained as previously described. 

The lift is determined solely from the pressure distribution on the cowl and 
vehicle undersurface assuming that the internal walls and side flow fences 
provide no contribution. Then, the lift is given by the integration 


L = 


(P-P ) Z(x) dx - 


(P-PJ Z(x) dx 


(18) 


1 


1 


The values of thrust and pitching moment contributed by the cowl and vehicle 
undersurface are - given by 


* Yv 2 ( % 

Th*= \ (P-PJ Z dy -. \ (P-PJ Z dy 


(19) 


* 

M ■ - 


2 ■ 


(P-PJ z (y-y) dy + \ (P-PJ Z (x-x) dx 


1 


1 


( 20 ) 


’ 2 (P-PJ l (y-y) dy - \ V 2 (P-P ) z (x-x) 


dx 


!\ 
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where (x,y) is the origin for the pitching moment. 


To assess the sidewall contributions to the thrust and pitching moment, 
momentum conservation integrals must be applied at the entrance and exit 
planes of the nozzle since the sidewall pressures are not explicitly cal- 
culated in the analysis, Before calculating these integrals, the exiting' 
and entering mass flux must be compared since slight discrepancies in mass 
flow can yield substantial errors in momentum' integrals. The mass flux is 
given by 

p v s i n ( © -0) . 



(21) 


where 6 g represents the inclination of the plane of integration with respect 
to the horizontal (i.e., for the exit plane depicted in Figure 1, 
e - tan * ((y -y )/(x u -x„ )). To perform an integration over the exit 

^ VpLn VqLq 

plane where properties are only known at the discrete points E -E.* the two 

U Q 

dimensional value of wave strength is linearly distributed between E points 
and properties are obtained at these interpolated points following the pro- 
cedures outlined in Section II. The properties are corrected to yield the 
correct existing mass flow as follows for frozen flow. 

Let R = represent the ratio of exiting to entering mass flux as determined 

by the integration of Equation (21). For frozen flows, the mass flux can be 
expressed solely in terms of the Mach number distribution 

( D(M)sin(0^e) 

• ?imf — Zrf y (22) 

'y 

where 

v 1 2 "(y + 1)/2(y-1) 

D(M) = M(1 + Y" M d ) (23) 

Then, at all points on the exit profile D (M) is corrected by division by R 
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D (M) = D(M)/R 


and the corrected Mach number is found by solving the transcendental equa- 
tion (Equation 23) for M by an iterative procedure; then, desired proper- 
ties are obtained by standard isentropic expansions 


P = P 


1 + ^ M 2 

1 + M' 2 


Y 

\ Y-l 


T 


= T 


1 + ^M 2 \ 

y ■** 1 ^ 2 

> 1 + 2^ M j 


For equilibrium calculations, at all points on the exit plane, the product 
(pu) is divided by R 


(24) 


(25) 


(pu) = 

I 

and the value of Av is sought that produces this change. Hence, an iterative 
procedure is employed wherein the value of Av at all points on the exit plane 
is perturbed to produce the corrected product (pu) employing Equations (7) - 
( 12 ). 


Then, having corrected the mass flow by modification of properties at the exit 
plane, the momentum flux is given by 



p\rsin(6 -e) cose + ( P-P ) sine 


[- 


sine. 


zdy 


(26) 


which for frozen flow, simplifies to 



yPM si n ( e -©) cose + (P-P ) sine,. 

[ S ; 2 i] Zd y 


sine. 


(27) 
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Then, the total thrust is given by 


Th 



( 28 ) 


where 1 and 2 represent integrations over the entrance and exit stations re- 
spectively, The thrust attributed to the sidewalls is then given by 




(29) 


The pitching moment is given by 
M - PM 2 - PM X 


(30) 


where 


and 


[f x (y-y) - fy(x-x)] zdy 
2 

pV sin(e -e) cose + (P-P )sine 

f = S s 

x sine„ 

s 

2 

. pV $1n(e -e) sine + (P-P )cose„ 
f = s s 

y sine 



(31) 


(32) 


( 33 ) 
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V. SAMPLE CALCULATIONS 

To assess the validity of the developed program, a comparison was made with 
the source flow characteristic program described in Reference (3). The noz- 
zle configuration and resultant wave field (as determined from the parametric 
design program (NOZD)) i s depicted in Figure (9). The data at various points 
in the flow field is given in Table I. A comparison of the pressure distribu- 
tions downstream of the cowl and undersurface expansion fans is depicted in 
Figure (10). The source flow characteristic program (NOZBOD) had the source 
located at x = - 20 and had the sharp turns at the cowl and undersurface occur 
over a distance Ax - .15, The source flow program calculates frozen chemistry 
via appropriate curve fits while the wave program calculates constant y frozen 
flows. It should be noted that the program NOZD required .4 seconds for this 
calculation as compared with 55 seconds for NOZBOD, indicating that the ratio 
of running times for frozen flow calculations is of the order of 140 to 1. The 
pressure distribution as calculated by NOZD shows excellent agreement with the 
second order source flow calculations of NOZBOD and the thrust values as cal- 
culated by the two programs agreed to within one percent. Additionally, de- 
picted in Figure (10) are the corresponding pressure distributions on the cowl 
and undersurface for the equivalent nozzle with no lateral expansion. The 
data at various points in the flow field for this calculation are given in 
Table II. 

To illustrate the parametric design process, a two dimensional calculation was 
performed for a nozzle subject to the following constraints: 


Entrance Height 

1.0 

Exit Height 

6.0 

Vehicle Length 

15.0 

Cowl Expansion 

10 


For these constraints, the cowl length was systematically varied from 3 to 9 
and for each cowl length the vehicle undersurface expansion was varied in one 
degree increments as depicted in Figure (3). 
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The initial conditions for this case are: 

P = 850 lb/ft 2 
M = 3 
Y = 1.4 
$ « 0 

and P w = 23.1 lb/ft 2 

The thrust is plotted versus vehicle turning angle with the cowl length as 
a parameter in Figure (11), and the lift is plotted in Figure (12). 



VI. 


CONCLUSIONS 


The numerical program developed should be a useful tool in rapidly assessing 
the affects of varying dominant parameters on scram jet exhaust nozzles, The 
results obtained by this program compare favorably with more precise analytical 
methods with a considerable savings in overall running time. While the current 
program cannot analyze all possible configurations, it is felt that once a user 
gains familiarity with the program, changes to accommodate other configurations 
can readi ly be made. 
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FIGURE 1. TYPICAL NOZZLE CONFIGURATION 
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FIGURE 2, GRID POINT CALCULATION 



FIGURE 3, VARIATION OF VEHICLE EXPANSION WAVE STRENGTH 



FIGURE A. F AND G POINT DOWNSTREAM OF EXIT PLANE 



' E„ 


■FIGURE 6. NO UNIFORM WAVE ZONE ON EXIT PLANE 
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FIGURE 7. WAVE DISTRIBUTION ON UNDERSURFACE 



Z 


FIGURE 8, WAVE DISTRIBUTION ON COWL 
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2. 1 3?0F>00 

2.1 8c0p-0 1 

1.1 066E*OO 

«1 ,745^F-0l 

3.717lr+00 

i.7*A?e-oi 


L c / 

2.9043E*00 

1„1R99F*00 

1 1452E+00 

3 # opB^f.oI 

4 # 2930F + 00 

3 # 02fl6E w 0 1 

' 

D 

6*8416E+00 

1 9 2nRflE*0Q 

1*34?1E+00 

1 * ?83a F»0 1 

e;,4797E + 00 

4 t 77?8E“01 


v3 

5 »49 1 BE* 0 0 

2.71 a!e *00 

1 .27A6E+00 

3 e 02 b<f-o1 

4»4l70F«00 

3*o2fl6E - *Ol 


C3 

4,l220E*00 

.7 2f.77E.0l 
€ 

1 206lE*00 
0 

.l # 745?F„0l 

3, 811 SF *00 

1.74F2E.01 




Mrt F POINT 
NO 0 POINT 


rv> 

Co 


i pT „ 

X 

Y 

z 

they a 

MACH 

ONijfO 

P 

T 

L EO 

5 . 0000 E +00 

-R e 8 l f 7 f« 0 1 

I . 2 F 00 E +00 

■• 1 , 745 '?F«- 0 l 

8 e 99 oRF +00 

2 . 1010 E -01 

2 6 1 B 47 F* 02 

6 o 7562 E+n 2 

F 1 

5 *S 189 E ^00 

• 5 * 7 o? 3 f -01 

1 . 2759 E *00 

“•1 » 1 P 9 ^E" , 0 1 

4 „l 2 l 3 F +00 

2 . 36 O 9 E -01 

1 * 81 30 F* 0 ? 

6 o 43 1 OF ♦ 0 ^ 

f 2 

8 f 3 * 7 lE + 0 ° 

1 12 a 7 E ^00 

1 4 l 74 E > 00 

1 283^^01 

R. 5195 F +00 

4 77 i 0 E„Ol 

3 4006 F 4 01 

3 9 R 67 E +02 

_.E 3 

X« 03 ^iE<fOl 

2 , 334 ^ 4.00 

ilsiaoE^oo 

l' 283 '.E-Ol 

B o 6 1 1 2 F + 0 0 

4 ’ 77 ^ 8 E -01 

3 I 0791 F 4-01 

3 l 875 lE *02 

p 4 

l*U 57 E 40 l 

2 . 8125 E +00 

i , 5 B 78 E + 0 0 

1 . 3 <? 9 qE -0 1 

B, 5691 E+Q 0 , 

4 .F. 574 E -0 1 

3, 222 If* 0 1 

3 , 9257 E ^02 

! E 5 

/ 

X* 5000 E* 0 l 

5 # llR 4 E *00 

1 7500 E *00 

1 . 9 n 5 iE. 0 l 

5 # 369 1 E + 00 

4 ,o 9 o 6 E „01 

4 . 0124 E +01 

4 # 1796 E ^02 


TABLE I DATA FOR THREE DIMENSIONAL SAMPLE CALCULATION 



xc 2 - 5 ,oo 0 of +00 

nN uv s 3 .n? 8 AF-nl 
nNUC s . 1 , 745 ?F- 0 l 
TAU = 1 1?R2**F-0l 

Ta u C = 2**f>laZ-02' 

TAURV a „I.09lRF.0l 

TAURC a 0 * 


; PT* 

l 

X 

Y 

7 

THETA- 

MACH 

' ONm?D ^ r 

I vl 

o. 

1 , OOOOr+OO 

1 .onnoE + oo 

3, 0?8aF-01 

4, 1372F + 00 

3* o2fl6E-0 1 

1 

1 ci 

0 • 

0 , 

l.onooE+OO 

-1 *745?F-0l 

8*6070F*00 

1.7452E-01 

A 

1.4264F+00 

5, OOO^F-Ol 

1 ,OOOOE*00 

0. 

3 ,0230r+00 

0 • 

8 

2 • 086 2 E ♦ 0 0 

2(,2?oOF"»0^ 

1 .OOOOE+OO 

- 1 .745^’F-O 1 

3 « 6070F + 0 0 

l,7^5 2 E-0) ; „ 

; c - 

2,81 3SE + 0 0 

] # ] 6t;4F+ 0 0 

l^OOOOE+OO 

3 # 0?8/, F«Ol 

4, 137?f*00 

3 # o 2 a6E- 0 1 

! n 

6*?253F+00 

uufyf.oo 

1 * OOOOE + OO 

1 ,283/,F-0l 

a • 068PF+00 

4.7738E-01 

j v3 

5«220*F>OG 

2.63]3f+ 00 

1 # oooof+oo 

3.02 Ra f-01 

4« 1 372f^ 00 

3,02 r6 E" 0 1 

L ~" c 3 ‘ 

3 • 97 03^ ♦ 0 0 

-7* Ono ? F“0 1 

iIooooe+oo 

-I .745^-01 

3.6070F+00 

I • 7^5 2 E-0 1 



.MO r POINT 




G 

i 

1 »o5l6E*ol 

?,6no 8 F*00 

l.OOOOE^OO 

1 •2934F-01 

k;*0 68 2E*00 

4.77 3 8E-01 

i 

i 

i.... 

. PT, 

X 

Y 

z 

THEtA 

MACH 

DNnpD P 


• ‘ EO 

5-0000^*00 

-8.8] c7f-0 I 

I * 

El 

5 « 6980F + 0 0 

•*4 , 6 ^ 7 8F«- 0 1 

1. 

i F2 

8 1527e + 00 

1 OIOOE^OO 

1 

! • E3 

l*f)773F+ol 

?;bp?5f+oo 

i: 

F^ 

lollROF+0 1 

?*8P63 f +00 

i. 

E5 

i,boooe + oi 

5 11 p4E . 00 

o T 

i 

* 


onooe+oo -i.745?f-oi ?.7908f+oo 

0000E+00 -9.00 5^-0 2 3*944] f +00 

onnoE+oo 1 283 af,o1 5 # o 6 8 2 f + oo 

nnOOE+OO 1 *2834F-0l F.068 ?f+o0 

nnoOF + OO 1 o 3 0 1 Qp-Ol 8*0797 f + 00 
.on ODE* 0 0 l # 937lE_0l 4 # 9023F + 00 


2*?1 rRE- 0 l 2,R?41F*02 7*2992 E*o2 

2 , 58992-01 2 .? 9 3 qf+ 0 2 6 * R 7 R 1 F ♦ 0 2 ! 

4 7738E_01 5 64302*01 4 # 6o74E + o2 

4 J 7738E-0 1 5* 6430F *01 4^60742*02 | 

^, 79 ? 2 e "01 5 *^ 685 f + 01 ^. 5 R 99 f*o 2 j 

4 # 50]6E w 0l 6^8504^^01 4 # 8698E*o2 .[ 




TABLE II DATA FOR TWO DIMENSIONAL SAMPLE CALCULATION 
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APPENDIX I 
PROGRAM DESCRIPTION 

A. Input 

Card 1 (Format 8E10.0) 

Column 1-10 P (initial pressure, Ib/ft ) 

11-20 Th (initial flow deflection angle, radians) 

21-30 EM (initial Mach number- frozen for frozen deck, equilibrium for 

31-40 T (initia^ r iempera^re, degrees Rankine) 

* For Frozen Deck 

41-50 G (specific heat ratio frozen) 

51-60 WGAS (molecular weight) 

61-70 PINF (external pressure P w , lb/ft^) 

* For Equilibrium Deck 

41-50 PHI (fuel-air equivalence ratio $) 

51-60 PINF (external pressure P , lb/ft^) 


Card 2 (Format 8E1Q.Q) 

Column 1-10 XV 1 (Xy of Figure 1, ft) 

11-20 YV1 (Y v ) 

21-30 XV2 (X v ) 

v 2 

31-40 XC1 (X r ) 

U 1 

41-50 YC1 (Y r ) 

4 

51-60 XC2 (X r ) 

4 

61-70 DYV (vehicle exit height Y y -Y r ) 

v 2 “2 

‘71-80 DNUC ( Avq-cowI turning angle-degrees) 
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Card 3 


Column 


1-5 (15) ICF-nurnfaer of different cowl lengths to be executed 

6-15 (E10.0) DXC-increment to be added to original cowl length X^ 


16-20 (15) 

21-30 (E10.0) 
31-35 (15) 

36-40 (15) 


OTF-number of vehicle turning angles to be run for 
each cowl length 

DTH-increment for vehicle turning angle in degrees 
IPR-print option 

IPR=0 prints properties on undersurface and cow! 

IPR=1 suppresses above printout 

IPGLY-polynomial print option 

IP0LY=0 prints polynomial coefficients of wall surfaces 
in form suitable for running program NOZBOD 
(Reference 3) 

1PGLY=1 suppresses above printout 


Card 4 (8E10.0) 


Column 1-10 

XSHF 

11-20 

YSHF 

Card 5 (8E10.0) 

Column 1-10 

AZ 

11-20 

BZ 

21-30 

CZ 


coefficients of lateral expansion quadratic 
Z(X)=AZ*X 2 +BZ*X+CZ . 

(where Z(X) is ratio of lateral nozzle extent 
at X station to that at initial station) 


B. Output 

The output is readily understood by referring to the figures in the report. 

If the F point or G point falls downstream of the nozzle exit station Eq~E: 5 » 
this is indicated in the printout by "No F Point 11 or “No G Point" respectively. 
If there is no uniform two dimensional wave zone* the exit line is not com- 
puted. The percentage thrust loss is based on the calculation of the ideal 
thrust using one dimensional considerations based on the entrance and exit 
area projections onto planes X=constant (i.e., ^ ex i t = (YV2-YC2)*(ZV2+ZC2) /2 ) . 
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PRECEDING PAGE BLANK NOT FILMED 

3. RHEQ-density fit in equilibrium deck p=RHEQ(P,T,<i>) 

4. ' FGAM-isentropic fit in equilibrium deck r=FGAM(P,T,$) 

5. FH-static enthalpy fit in equilibrium deck h-FH(P,T>$>) 



APPENDIX II 


LISTING OF FROZEN FLOW PROGRAM 
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PROGRAM NOZO (INPUT .OUTPUT •TAPEK«lNPUT t TAPE6aOUTPUT) 

COMMOm/A/PI ( 3) ,Twi (3) * E M 1 (3) tAl (3) *G1 (3) 

COMMOM/SHF/ xSMF.ySHF 
rOMMOM/POl./ I PR * T POLY 
cOmmom/thRmax/thrmax 
■ n I MENTION T 1 { 3 ) 

COMMOM/B/XA,YA,XB f Y 8 ,XC, YC f XD»vD» 
lTHA*TH8*THC9THOtEMA?EMB*EMCf EMn,xMUAf xMuBtx^UCfXMUD 
COMMON/C/ONUVtONUCtXCl «YCl *XVl *YV1*DNUL 
c 0MM0Kf/ n 3/ AZ,rZ* c Z 
COMMOM /V/ AVfBV«CV 
OOMMOm/P/ PINF 

DIMENSION P2(I0) «TM2(10) »EM2(ln) 9T2(1 o) *62(10) *X2(10) «Y2U0) 

l*Z2O0> 

DIMENSION HOL(U) *H0LE<6) *DNuE(10) 

DATA hOLE/2HEo* 2HE1 *2HE2 # 2HE3t?HE4*2HF5/ 

DATA h°L/2HV19RHc 1*2HA *2Hs i2H0 *2HV3 , 2HC3* 2HF * 2HG ?2HV4/ 

GETz <yS) =AZ*x5*XS*BZ*x5*CZ 

98 FORMAT (7E13. 5) 

99 FORMAT (/) 

100 FORMAT (AFiOpO) 

PEAO(*i*i 00) PI (I) «THI (1) 9 EM1 (1) ,T1 (1) « 6 I (1) ;wGAS*PINF 
DO 106 1=2*3 

Pl(I>aPl(l) S TH1 ( X ) =TH1 (!) % fMI ( I ) »F«1 < 1 ) $ G 1 < I ) =G1 ( 1 ) , 

Tl(I)sna) 

106 CONTINUE 

DUM =5 . + (G1 < 1 > -] * > / 2 • ^ F. M 1 { 1 } 

PTOTkpI ( 1 ) ( G 1 <l)/<Gl(l>-t*>) 

TTOT = Tl ( 1 ) *PUM 
RC, = 49aOO. • 

RGASbrG/WGAS 

TAUPnaO, 

READIb^ioO) XVl^YVl »xv2’ XC19VC1 5 XC2?DVV?0NUC 
PEAD(B9B921) IcFoDXC* JTF iDTH 9 tPR » IPOL y - 

5921 FORMAT (r59ElQ*0*I5«ElO.O f 2l5) 

PEAD(BtinO) XShF,YSHF 
ONUCedNUC/57,3 

RFAD (5*1 00) AZ 9 BZ ? CZ 
ZCl5AZ"XCl tt XCl^B7^XCUCZ 
ZVi=A7»XVl*XVl+BZ*XVl*CZ 
ZC2=GF.TZ (XC2) 

WRITF(6 t 5922) 

5922 FORMAT ( 1 W 1 /////) 

WRITE (6^1020) 

1020 FORMAT* 2 0 XttF R 0 Z E N FLOW N 0 7 7 t E r> E S I G N«) 
WRITE (6*400 ) P'l f 1 ) * THl ( 1 ) * EMI f 1 ) 9 T 1 ( 1 ) 9 G 1 ( <] ) 9 WGAS*PImF 
400 pQRMAT * ///BOX, ^INITIAL PROfr ILF*//3 X*PResSURE#7X#ThET A*8X*MACH* 

16X*TemPERATURE*5X*GAMma* 7X*M0L WT#8 X #e> I NF#/7E 1 3 , 5 } 

WRITE (6*5930 ) X V 1 9 Y V 1 9 X Vg « XC 1 * vC 1 9 XC2 9 DY V 
59 30 FORMAT {///6x*xvl# l OX*YVl*lOX*Xv2*lOX*xCl*lOx*YCl*lOX*XC2*lOX*DYv*/ 
I7El3 # q///> 

WRITE(6,5923) xSHF.YSHF 

5923 FORMAT <9X*M0MENT aXIS*/12X*X =#£ 1 3 . 5/ 1 2X*Y s*e13,5///> 

WRITE ( 69 1001) AZ*BZ«CZ 

1 00 1 FORMAT (9X*LATERAL EXPANSION FQiiaTI0N#/9X92?HZ(X) S aZ#X##2*BZ #X + CZ 
‘ X/12X«aZ ■*EI3.5/T2X*BZ =#E13.5/12x*CZ =#Ei3.5) 
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XC22=yC? 

THC1=tH1 (3>-DNmC 

DO 5 00 0 I XC = 1 * T CF ... 

XC2=XC22 + FL0AT flxC-lj^DxC 

CALL THM (EMI < 1 ) tPl < 1 ) «G1 (1 ) pPIMF*YV 1 »YC1 •DYVtXV3f THRMAXf XC?) 
WRITE <6,6364) JHPMAX 
6364 FORMAT (lHl3iX^TDEAL THRUST =tfET2 t 4) 

YC2=YC1+TANCTHC1 ) # (XC 2 -XCU 
YV2=Y^2+qY V 
DO 6 0 0 0 jTsltJTF 
THX=1.E+I0 

THJT sFLOaT ( JT-1 ) *DTH 

DNUV=aTan( (YV 2-YVl)/(XV2-XVn )-THl (3) JH JT/57.3 
CALL dM(FM 1 ( 1) t OMUC,EMCl,Gi ( n ’) 

CALL pM(EMl (1) ,DN«uV»EMV1*G1 ( 1 ) ) 

IT2 = 1 

1 0 00 CONTINUE 
CALL CHAR 
DNUAsn* 

nNUD«nNUV+QNUc 

7A=GEtH{XA) 

7P = GftZ ( X0 ) 

7C :sG ETZ(Xc> 

7D=GfjZ (XQ) ' 

pat=o, 

THC3c T HC! . - 

CALL PM(£TM1 (1) ,0NUC»EMC3,G1 (\) ) 

CALL FIX ( XC 1 *V C l*THcl«0.,X 9 ,YB9THB*XMU8»Xc3,Yc3*THc3*xMUc3feMc3« 
101 ( 1 )) 

7C3=GfTZ (XC3) 

IF(Xc3<»LT.XcH) GO TO I 77 
TAuCsO, 

XC3=XC2 

DNUEOsDNUC 

THC2cTHC1 

THEOsriHCF 

CALL PM(EM1 ( 1 ) t DNUE 0 , EME 0 « G 1 (1) ) 

CALL FM3D(EME0,EME0N,XC2 9 Gl ( 1 ) ) 

EMEQsfrMEnN 

XMUEOsASTN ( 1 */FMEO) 

GO IQ 175 
177 CONTINUE 
IT9 = l 
RATs.1 

176 CONTINUE 

XMUGOkXmuB+RAT^ (XMUO-XMUb) 

THBD=THb+PAT* (THn^THB) 

XT=Xr*R A T# (XD-X8) 

YTsYr+RaT* ( Yq-yB 1 
TaUc=RAT*D n UV 
TE(Rat«GT.IoO) <30 TO 6000 
DNUW-nNuc ♦ T AUC 

CALL pM(EMl ( 1 ) ,ONUW»EMW,Gl (1 ) ) 

THW=ThC3*TAUc 

CALL Fix (XC3» YC3*THC3*XM(jC3*xT *yT , THBn 9 XMUBD 9 X'W 9 YW » THw , XMIJW 5 EMw , 
1 G 1 ( 1 ) ) 
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DNUEnsDM uc*2# <*yauc 

THc2=tHc1 

THEO=tHC? 

CALL pM(EMl (1) ♦ DNUE 0 t EM£q « G l (1) ) 

CALL rM3n (EMEO,EME;ONf XC2*G1 (1 ) ) 

EMEOsfMEoN 

XMUEoskASIN < 1 J 

$LWC?= (TAN (THW-*XM|jW) + TAN ( THEO-YMUEO ) ) /2. 

" yC2T B yW*SLWC2# (XC2-XW) 

FP9=Yc2“YC2T 

TF(ARS { ER9> oLT.l.E-03) GO TO I 75 

CALL ERROR (4*lT9,RAT,ER9 # l.li»PATl»ER9l) 

GO TO 176 
175 CONTINUE 

THV3*stH 1 ( 1 ) ♦ DNlJV 

CALL pM(EMl (1) t DNuVtEMv3,6in) > 

THV1=tHV3 

CALL FIX (XC#Yc*THC»XMUCtXVltYVl f THVltO*#XV3*YV3»THV3«XMUV3*EMV3 f 

l a l ( 1 ) ) 

A V a Y V 3 S 8V=TAM(THV3> 

CV= ( YV2-AV”BV^ (XVP-XV3) ) / ( X V?- X V3 ) **£ 

7V3=Gf t 7 (XV3) 

THV2 = /\TAN oCV# (XV2-XV3) ) 

IF(Thv2pLT«0*> GO TO 462 

TAURV=THV2«THV3 

IT?=1 

TTAUs:^ 

P A T = 1 , 

545 CONTINUE 
530 CONTINUE 

xz=Xc< R at* (Xd-xC) 

YZ=Yc+RAT* ( yd~yc> 

TA(J s R AT <>f)NUC 
THZ=ThC+PAT* (ThO-THC) 

EM7=EMC+RAT# (emd-ENC ) 

X.MUZsflSlN ( 1 #/£mZ ) 

T.Hy = THV3wTAU 
DNUYsh^UV+TaU 

CALL pM(FMl (1) fOMUY*EMY»Gl (1 ) ) 

CALL Fix (XZf YZ,TB7»xWuZ>xV3oYV3,THv3i X M(jv3*xYtYY»THy*x M UY*E^Y9Gl < 1 

m 

DNUE5=ONUV*2.*TAU*TAURV 

THE5=tHv? 

CALL PM(EM1 ( 1 ) f 0MUE5>EME5fGl (1) ) 

CAll FM30 (EME5,EmE.5Ni XV2»G1 (1 ) \ 

XMUE5aASIN'(l«/EMF5> 

SL2Ts(Tan (ThY+XMUY) ♦Tan (THE 5+XMUE5) ) /?t 
YV7>T = y Y*SL2T# { XVP-XY) 

IF ( IT7*EG» 1 •ANn. YV2T.GE* YV?) Go TO 54? 

IF(APS(E7) .LT.UF«03)GO TO 540 

IF (OMuCeEO* u * ) GO TO 540 

CALL fRROR (5*ItZ*RAT.EZ» „ 1 uRATl .EZl ) 

GO TO 530 
542 iTAUof 
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540 CONTINUE 

ONUEjsONUV ♦DNUC 

CALL pM{fMl (1) f DNUE2*EME2fGl (In 

THEH=THi (1) ♦DNNV^ONUC 
THF = TnE?> 
p-Mr = EME2 

CALL Fix (XC3*YC3*THC3*XMuC3*XD,YD»TH0,XMUD f xF,YF,THr»xMUF»EMFt 
1G1 (1 ) > 

SLE K (YV2*YC2) / (XV?~Xc2> 

THE-ATAN { SLE) 

EMFTs.pME2 

CALL rlX (XC3* YC3*THC3«XMUC3*XC?»YC2*THE?0» f XFT*VFT?TNFiXMUFT»EMFTo 
lGl (1)) 

IF (XF,GT,xPT) GO TO 300 

CALL FIX (XF*YF,TMF*XMUF9XC29YC?*THE90, t XE2 9 vE2 9THE2 * XMUE2 * E^ES t 
lfil (1)5 

ZFsGEtZ(XF) 

SLF = taN ( tHF-xM||'F ) 

CALL GEM<XF*YF*SLF»XC2tYC2*SLE.XEltYEl) 

TTEl=0 

36B CONTINUE 

PAT= <yE]«YC 2) / (YE2-YC2) 
nNUElaf)NUV*ONUC+ (1 *-RaT) *TAUC 
CALL pM(CMl(l) t DNuEltEMEi,Gl(rn 
THF 1 JsyN 1 ( 1 ) ♦ DNl.iV-DNUC- ( 1,-RAT) *TAUC 

CALL fIX (X c2 ? Yc 2 9TH£»0. * XF * Yp • THF ♦ XMU f » X£ l * Ye 1 » TheI * XmUfI »fMe1 * 
lGi(lp 

IF(lTFl,EO«l) go to 359 
ITEl=l 
GO TO 368 
369 CONTINUE 
GO To 3 in 
300 continue 

CALL FIX (XC2»YC2tTHEt 0 . tX0fY0*THD*XMUn*XE2»YE2*THE2«XMUE2tFM£2t 
loidu 

SLCE-TAN ( THC3+XMUC3 ) 

CALL f?EM(XC3*Yc3fSLCE*XC2 9YC2 9FLE»XEl*YEl ) 

I TE 1 ~0 

38 l CONTINUE 

RAT=(YE1-YC2)/{YE2-YC2) 

DNUFUDNUOTAUC + RAT* (DNUV-TAUC) 

CALL P M ( F M 1 ( 1 ) t DNuEl,EMEi,Gl (1> ) 

THF1=tH1 (1 ) -DNUOTAUC+RAT* (DNUV/-TAUC) 

* CALL Fix (XC3« Yc3fTHC3»XMUC39XC2t Yc29TnE90 o ,XE1 *Ye1 *THfI 9XMUE1 eFMFl 
1 t G l ( 1 ) ) 

TF(lTFl.EQ.l) GO TO 380 
T T E 1 ” I 
GO TO 381 
380 CONTINUE 
310 CONTtmUF 
THE3-THEP 
Dmi |E 3sDwi (E 2 

CALL PM (EMI (1 } ,DNUE3,EME3*G1 (lj) 

XMUE 3 = A 5 1 N ( 1,/fMf3) 

TF(ItaU.NE,1) GO TO 311 

XG~X Y $ YG--YY % YHG-THY S EMG=FMY $ XmUG»XMijY 
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xgt=xg 

GO TO 3?08 

311 CONTINUE ‘ 

THG^TmE? $ EMG=EME3 

CALL rlx {XO*YD, THDfXMUD»xV3, yV3,THV3 t yMUV3 f xG,VG,THG« xMUG,EMG, 

161 (15) 
pMGT =EMp3 

CALL Fix {XC2 t YC2»THE,o. *xV 3 #YV 3 .THV 3 *xMUV 3 , xGT»Y 6 T,THG,XMUGT»EMGT, 
161(D) 

TF(XG,GT.XC5T) GO TO 350 
3588 cONTtmUf 

CALL FIX (XC2» YC2*THE*0* *XGfYG»tHG»XMUG*XE3»YE3»THE3«Xv|UE3«EM£3 9 

loi < i n 

7G=GETZ(XG) 

IF(YE3*LE.YE2) GO TO 359 ‘ 

SLE /+ =TAN ( THG + XmUG ) 

CALL GEM{XG»YG,5LE^»Xc2*YC2f SL f»XE^*YF 4 ) 
tTa = o 

358 CONTINUE 
RAT=(yE 4-YV2)/ (YE3-YV2) 
nNUE*=DNUV^[)NUC + (1 *-RaTMMTaURv + TaU) 

CALL pM(FMl (1) ,DNUE4,EME4*G1 (In 
THEAerHl (1 ) ♦ONIJV-ONUC* U . -RAT ><M TaURV + TaU) 

CALL Fix (XGp YG*THG*XMU6tXC2§YCpf THE *()• 9 XE A « YE4 1 THE A * XMUE4 • EME4 ♦ 
161(1)) 

IF ( I T 4 *E0* 1 ) GO TO 359 
T T 4 = 1 
GO 358 

359 CONTINUE 

TF ( I TaU, EQ# 0 ) GO TO 360 
SLV32- ( YV2^YV3) / ( X V2-XV3 ) 

SLGsTaN(THG*XMDG) 

CALL GEM(XV3fYV3*SLV32»XG»YG#SLG*XV4tYV4> 

TTV4 = i 

362 yV4 = AV + BVtt (XV4p,XV3) + CV#(XV4«,XV3)**2 
THV4 «aTAN (BV+2,*cV* (XV4-XV3) ) 

TAURV4 b THV4«-TH\/3 

t)NUV4=0NUV*f .*TAU^TAURV4 , 

CALL pM (pMl { 1 ) .L1NUV4 ♦ fMV4 1 Gl < 1 ) ) 

CALL f m 3D (EMV4 t EMV4N, XV4 f G{ ( 1 ) ) 

FMV4«fMV4N 
XHUV4-ASTN < I./fHVA) 

YV4T=syG*,5o ( TAM<ThG + XMUG) ♦ T AN ( THV4 ♦ XMlJ V4 ) ) # (XV*”XG) 

FRV4 B ‘ ( YV4T-YV4)/YV2 

IF ( A p 5 ( E R V 4 ) • L T • 1 • E * 0 3 ) GO TO 361 

CALL fRROR <36l , iTV^iXV* *ERV*» 1. 1 >XV 4 1 «ERV*l ) 

7 v 4sGfTZ ( XV 4 ) 

IFUTv4.EQ*l ) GO TO 361 
GO TO 36? 

361 DNUE4aDNUV4 

TF(YE3*LE* YE? ) GO TO 460 

CALL pM ( FMl { 1 ) ,DN!JE^ 9 EME 4 »G 1 (In 

THE4 = tHV4 

call Fix <XC 2 tYC 2 *THEfO. *XV4*YV4fTHV4?XMUV49XF4»YE4«THF4tXMUE4* 

. lFNF4,r,l (1 ) ) 

GO Tq 360 
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350 CONTINUE 

1 T&k o 

CALL FIX (XO>YD,THO*XMUD*XC2*YC?,THEtO ft tXE3 * yE3 * THE3* XmUE3 *EME3« 
101 ( 1 )) 

IF<YEVLE.YE2) go to 460 
SLE4»TAN (THV3-XMIJV3 ) 

CALL gEM (XV3,yv3,sle4,Xc2*Yc2*FLE*XE4 9 YeA) 

371 CONTINUE 

rAT« (YE4-YV2) / (YE3-YV2) 

DNIJEAsONUV + TAU^RAT 4 * (DNUC-TAU) 

CALL PM (EMI (3, ) % DNUEA,EME4iG1 1 1 \ ) 

THE4stHi ( 1 > + DN!|V**TAU«RAT# (nNUC«TAU> 

CALL pIX fXC2?Yc2*THE90,#XV3fYV3*THV3»xMuV3tXE4»YF4,TH f r4tXMUE4fEMF4 
1 * 61 ( 1 )) 

TF(IT4*E0.1)G0 TO 370 

TTA = 1 

GO To 371 
370 COMTTmUE 
360 CONTINUE 

TF (YE3*LE. YE2) GO TO 46n 

X ? ( 1 ) r= X c 2 $ X2 ( 2 ) s=XE 1 S X2(3)ey E 2 $ X2<4 )=Xf3 $ X2(5)sXEA 

X 2 ( 6 ) = X V? 

y2 ( 1 ) s YC2 $ Y2(2)3Y£1 $ y2(3)=yE2 5 Y?(4)sYF3 $ y2(5>-yE4 
Y2(6>-Yv? $ EM2(1)sEME 0 $ EM2(?>=EMEl * EM2^3 )=EmE2 $ £H2(4)=EME3 
EM2(5)=EHE4 $ fM?(6>=EME5 $ THr(1)=THfO $ Th?<2)sTHf1 
TH?(3)=THE2 $ TH?(4)=THE3 * TH?(5)=THf4 $ Th2(6>=THE5 
ONUE ( i ) sONUEo * DNUE (gJaONUEi < DNUE ( 3 ) sQNUF ? $ DNUE ( 4 ) =DNUE 3 
pNUE < 5 > shNU^A * nNUE (6 ) = d NUf 5 
INDEXsO 
1 BO CONTINUE 

00 Ifio 7 = 1*6 

nUM- 1 # + (cl 1 1) *1 « ) /2 « *£M2 ( I ) -»>*2 

02 ( I ) -PTOT/DUM*},# (G1 <1)/(G1 ( 1 ) -Y * ) ) 

T?(I)sTTOT/DUM 

72m=AZ tt X2(I)**2*B2*X2(l) *cZ 
160 CONTINUE 

THSlsi .5707963 

IF (XV) .NF.Xcl) THS1 s ATaN( ( YV 1-YC* > / < XV UXci) ) 

THS2=aTAM ( (YV2-Y02) / (XV2^XC2) ) 

7I = .5*(GETZ(XCl ) ♦GETZ (XV] ) ) J . 

XMASS 1 S P 1 (1) *EMl ( 1>/5QRT(T1 : 1 ) ) *SIN(ThS 1-TH1 ( 1) )/SlN(THSl)*Zl* 

1 (YVUvCl ) 

OUM = Gi (n^EMi (1 )#*2#SIN(THS1«Tni <1> )/SIN(Ths1) 

U s 1 c-plNF/Pl ( 1 ) 

THXl=pl ( 1 ) *<U*DUM*COS (TH1 U ) ) ) *?I # { Y V 1 -YC 1 ) 

CALL xTH (X2»Y2,TH2 f UNuE,EHl 9 G1 ,PTOT, TtOt» TH fi2 9 THxi * xMaSS1 » THX) 

ISHsO 
GO TO 461 

460 CONTINUE 
ISH*j 

WPITE(6»99ll 
99 I fORM a t(* Ye3.Lf* 

46 1 CONTINUE 

WRITE (6* ] 004) XCp.DNUV tONUC? TAtU TAljC 9 T aURV 9 TAUPC 

1 0 p4 f ORMat< 32X*X C 2 s*p;12.4/32X*dNUv 1 2 , 4/32X#pNUc = *El2«>4/ 

]32X#TaU s*E 12,4/ 32X&TAUC = *E12 9 4/32X*TAlJRV a#El2,4/32X*TAURC s 
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l#El2*4///> 

WRiTp f 6 f 5003) 

10 03 FORMAT < 1 X* p T,* 6 X*X*llX*Y*llX*Z* 8 x*THRTA*7x*MACH*8x*Dnu2D*/> 

WRITE ( 6 * 3 0 02) HOL { 1 ) * X V 1 * Y V 1 ♦ 7 V \ v THV 1 1 fMV 1 9 DNUV * HOL < P ) *XC1«YC1 • 

17Cl *ThC 1 tEMcl , nNi)C*HOt (3) f XA*Ya ,ZA»THA*EMA.nMUAfHOL (4) ,XB^YRiZB 9 
lTHBtEMB*DNUC*MnL (5) * XC t YC * ZC * TwC *EMCf DNUV, HOL ( 6 ) * XD , Yn * ZD * THD 9 £MD « 
1 DNUD»hOl (7) « XV 39 YV 3 >ZV 3 *THV 3 #EHV 3 t DNUV* HOL (a) ♦ XC 3 ♦ YC 3 * ZC 3 ? THC 3 » 
lEMC3*nNur. 

1002 FORMAT ( 1 X A2 9 6Ei 2 5 4 ) 

IF (Xp.GT.XFT) WRITE (6tl005> 

1005 FORMAT <20X*NO F POINT*) 

IF(XF.LF # XFT) wRTTE ( 6 f 1003) HOL (9) , XF , yF f ZF 9 f MF , EMF 9 ONuD 
TF (XG. 6 T # X G T ) WRITE ( 6 * 1006) 

1006 FORMAT <20X*NO 0 POINT*) 

IF (X6.LE.XGT) WRITE <6 , 1 002 ) HOL ( 10 ) , XG, YG , ZG f THG f EMG,D nUD 

IF ( I TaU 0 EO« 1 ) WRITE (69 10 02) HOL < \ \ ) * X V4 $ Y V4 * ZV4 » THV4 9 EMv/4 » DNUV4 

WRiTf( 6 *1013> 

1013 FORMAT ( ///) 

TFUSH.EO.D GO TO 7384 
wRlTE(6tl00M 

10 07 FORMAT (1 x^PT.^fcX^X* 1 ! X^Y^l l X X*THEy A«7 x*mACH*H x*DN| j2U* 1 0 lX 

WRTTf (8 f 1010 ) (hole a > »X2 (I> *y2(I) *Z2(i> *Th 2(I) *EM2(I> tONUE(I> * 
1P2(I) f T2{I) *Isl tfe) 

1010 FORMAT ( 1 XA2t8El 2.4) 

7384 CONTINUE 
XI=XVM 
THI-THl <\) 

CALL pO( Y{XVlvYVltTHI^XV3*THV3,XV2«TH\/2*XCl*YCl»XC2tTMC2«XT) 

CALL THRUST (XVl^YVitEMVl • XV 3 * YV 3 f THV 3 , XV 2 » XV 4 « TaU*DNUV» 

lXCl«Ycl*EMCl »XC3*YC3*XC2*VC2» T aIJC *DNUc ? PTOT * ITALUTHX) 

WRITE (8,99) 

WRITE-(6,99) 

OO To 6000 

462 CONTINUE 

WRTTF (6*992 ) 

992 fqRMAt (*THV2«lt*0*) 

450 CONTimUe 
6000 CONTINUE 
5000 continue 
END 


\ 
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SUBROUTINE CHAp 

CQMMOnj/A/PI (3) *THi (3) ? EMI <3> »A} (3) *61 (3) 
CO^^QM/B/xA^YAtxBfyBfxC^Y^’X^’V 0* 

1 THA f THb , THC , THO , EMA * EMB t EMC , EMn , XMUA , XMUB , XMUC * XMUD 
COMMOM/c/l)HUVf ONUC * XC 1 * YC1 » XVl • YVj ♦ UiMl.lL 
G*G l ( 5 ] 

x mu l = a s i w ( 1 ,/emI r l ) ) 

E M A = E Mi (] ) 

. THAsTHKI) 

CALL FIX fXCl t Ycl *THl ( 1 ) * XMl.ll tXvl * YVl * THl (1 ) , XMU1 * 

1 X A 9 Y A • JHA » xMyA »£ma »G) 

THR= THA-OWUC 

CaLL pMOrM 1 (1 ) ,uMUC’E m B*6) 

XMUfcj] ,=AS tN (1 # /p.M81 ) 

CALL FIX f X C i , ycl f THd * XMUtfl * XA * YA * THA t XMUA t XG * Y0 *THB* XmU 8 tEM8 *6) 

THC = 1>U ( 1 ) ♦D'^UV 

CALL PM (EMI ( 3 > ,DMUV*EMC*(5) 

pHC 1 = r M C 

XMUCl-ASIN ( 1 ,/FHC] ) 

CA LL PIX (Xa» Ya,thA*X h uA?XV 1* YVl * THC 9 X M UC 1 9 XC 9 YC » jHC » X M yC » E M C * 6 ) 
60 D N UP=D N UV*U N UC 

CALL PH (KM 1 ( i ) .DMUD»EMD*G) 

Thd = T h 1 (1 ) + N U V '* p MUq 

CALL F 1 x ( XB * yU * THb « X M UB ♦x ( - : *YC*THC 9 x !V, UC?x^*YP*^^^'? | XMUD,FHD*6) 

RETURm 

END 
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SUBROUTINE THruST (XV 1 * YV 1 *EMVl • XV3 1 YV3 * THV3 »X V2 * XV4 * TAU *DNUV • 

1 XC1 *Ycl ?FMCI *XC3»YC3*XC2»YC2 *TaUC*DNUC»PT 0T, ITAU,THx) 
COMMOm/a/PI (3)fTHl<3)tEMH3)tA1<3)»Gl(3) 

C0MMON/D3/ AZtftZfCZ 
COMMON /V/ AV*RV,CV 
COMMon/P/ PINF 
rOMMOM/SHF/ XSHFeVSHF 
COMMOm/POL/ I PR f i poly 
COMMOn/ThRMAX/THRMAX 
COMMOn^PmPM/ PmOm 
100 FORMAT ( BF 1 3 • 5 ) 

GFTZ(Y5)=A2' i!, X50Xc;4.BZ tt X5 + c2 

PI (EM5) sPTOT/U .♦ (Qi 1 1 )-l.) /2 *#eM5**2)* # iQUl) / IqHV -U) ) 
XLTF=n« 

X T H = 0 • 

XMOM=o • 

Dx=(Xv 2 -XVl)/50 e 
N= (XVl^XVl ) /OK* 1 • 

OXN* (XV3-XV1) /FLOAT <N) 

N = N+1 

CALL PM(FM1 (1) ,pfsl U v*EMx!*Gl ( 1 )> 

XA = XvT % Y A = Y V 1 % ZA=GETZ<XA> 

CALL EM3D(EMXl f EMAfXA*Gl (1) ) 

PAsP? (EMA) 

PX-PA ♦PlMF 
PX=PX/PTNF 

IF (lPp*£Cwl) GO TO 50 
WRITE (6, 110} 

1 JO FORMAT i///35X*\/EHXCLE UNdE PSURf ACE*/7 1 2 X# Y* 1 ?_ X*Z<* qX*P/P INF* 

1rX#Ma c H#qX#LIFT#8X#THRUST # 7X<*MomENT#> 

WRITE (6, 100) XA»YA*ZA f PXfEMAfXLlFf XtH f XMOM 

50 CONTINUE 

no 10 I=P«N 
XB=XA*OXN 

Yr«Yvi * (XB-XV l) /<XV3*XV1 ) #<YV3„YV1) 

7B=GFT2(yB) 

CALL E M 3D(EMXl,EMB*XB f Gl (1) ) 

Pr=P7 (EMr) 

nUM* (PA+PB } (Zn*7 A ) /4 # 

OXLFs dumMXB-XA) 

DXTH-oUM^ABS ( YR mYA ) 

YBB= ( YA+YR ) /2o-YSHF 
XBB=- ( (Xa+Xp) /P.-XSHF) 

0'xMOM s ~YBB»DxThUxBBttDxlF 

XLIF = xLTF' + OXLF 

XTH=XTH*nXTH 

X MOM sy MOM* OX MON 

X A = X R $ YA = YB 5; ZA = 2e $ PA = PR 

EMA^emH 

PX=PA+PTMF 

PX bPX/PINF 

IF ( IPR.EO. 1 ) GO TO S’ 1 

WRITEfGtl-OO) Xa*Ya*Za*PX,EMa«XlIF*XTH,XM0M , 

51 CONTINUE 
.10 CONTINUE 

N= (XV2"XV3> /OK* 1 • 
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DXNo(XV2-XV3)/FL0ATCN) 

NsN^ X 

DO 20 I=?fN 

xb=xa+dxm 

YRbAV*«V* (XB-XV3) ♦ CV 4 MXB-Xv3)*»2 
7 R= 6 PTZ (XB> 

THR=AtAN (BV+ 2 .#CV* (XB-XV 3 ) ) 

TAUR8-THFuThV3 

pATBr; (XR-XV3}/{XV3-XV3) 

TF(lTAU # pQol) p A T 8 55 ( X B “ X V 3 ) / ( X V A • X V 3 ) 

TP (RAT 8 .GT. 1 , ) RAT 8 = U 
r)NUe = nNUV*TAURR + ?,*RATB*TAU 
Xl=Xv3 S X2=XV? 

IF(ITaU.EQ* 1> x 2 = x V ^ 

RATAs ( (Xl-XCl ) / (XB“XCl ) -1 0 ) / ( (Xl-XCl ) / (X2*XC1 • ) 
IF (RatA,GE.1.) RATA-1. 

T AI.J X = A T A N { R A T A # T A M ( T A U ) ) 

DNUB=D^UV*TAURr«?,^TAUX 

CALL pM(pMl (1) f oNUB»EMX»el (1) ) 

CALL F m 3D(EMx *EMB* X 6 *G1 ( 1 ) ) 
pRrPZ (EMR ) 

DUMs (pA+PB) « (2R+ZA) /4* 

nXLF= oum*(Xb-xa> 

OXTHooUMttABS (Y r-Ya) 

YRB= (YA*YB> / 2 *-VSHF 
y R B a «• ( ( X A + xB)/?,-x 5HF > 

QXM0m = -YB8* > DXTH + XB8#DXLF 

XL IF=xLTF+DXLF 

XTH=^TH+nXTH 

XM0M=xM0M*DXMOM 

Xa=Xr $ YA"YB % ZA= z e £ PAsPr 

FMA=fmB 

PXsPa*.PINF 

PX=py/PiwF 

TF(IPP*FQ. 1 ) GO TO 52 

WRITE ( 6f 1 00 ) xA*yA*zA*Px?E m AjX|,I.F?XTM*xMOM 
52 CONTINUE 
PO continue 
XLVsXLlF 
XTHVbxTH 
XM0MV=XMOM 
WRITE ( 6 * 101 ) 

101 FORNAT(Z) 

X A = X C 1 S Y A = Y C I 5 ZA=GETZ<XA) 

XLIF-Ob 
X.TH = n . 

xmom=o* 

Nr { XC?- X C1 ) /Dx + 1 • 

DXN-(XC3-XCl)/FL0AT(N) 

N=N*} 

CALL pMlEMl (1) *DNijCiEM Y 1 ( 1 ) > 

CALL FM3D(EMYl # EMAtXA,6l (I) ) 

PA-FZ(EMa) 
p X r P A + P I ME 
pXrPX/PTMF 

IP ( I pp » pQ 8 1 } go TO 53 


{ 



TR 188 


Page A2-11 


WRTTFT/6*lll) 

i t i FORMAT <///50X*COWL# /7X#X*l2K#Y#l2X*Z*pX tt P/PlNF# 

1pX^Mac H ^9 x *L I PT^bX^THRUST«7X*MoMENT^) 

WRITE? 6, 100) xA»YAf ZA*P x*EMA,XI IF,XTH,XM0M 

53 CONTINUE 

DO 3o I=2*N 
XR=Xfl+DXM 

Yb s ^C 1 4 (XB-Xcl) / (XC3-XC1 ) # (YC3-YC1 > 

7B=0FTZ ( XB> 

CALL pM3 D (EMYltE:MB*XB*Gl <1> ) 

PB=P? (EMB) ' 

DUM- (pA + pB) «- (ZR«-ZA) /4, 

DXLF*s«DUM*(Xb-xa) 

DXTHsnUM*ABS (Y r-YA) 

" ' ' YBB S <YA + Y0} /2,-YSHF 

( (XA+XH) /?»«xSHF) 

DXMOm = w YBB^DXTH* XBS^DXLF 

XLIF»XLTF*DXLF 

XTH»XT H +nXTH 

XMOMsxMOM+DXMOM 

XA = Xp $ YA”YB % 7A”2B $ PA=PB 

PXcPa+PINF 

PX=PX/PTNF 

TF(lPP*FO*l) GO TO 54 

WRITE <6, j 00) XA» YAt 2 A*Px*EMA,XI. lF,XTH,XMOM 

54 continue 
30 CONTINUE 

M» ( Xcp*XC3 ) /DX + 1 « 

IF < N*pQ« 1 ) GO TO 60 
DXN=<xC2„xC3) /FLOAT ( N) 

N*N+1 

pO 40 I»2»N 

y9sxA*O x M 

YB«YC3* <XB-XC3)/<XC2-XC3)*(YC2«YC3) 

7b=6FtZ(XB) 

PaT= (XB-XC 3 ) / <XC?^ X C3) 

DNUB,oNUc+RAT* ?# *TAUC 

PATAs ( (XC3”XV \ )/ (X B"*XVl ) ~1 9 5 / < ( XC3"XVf > / ( X C2 W XV 1 > ~ 1 e > 

TAUYsaTaM<RBTa*TaN<TaUC> ) 

pNUB=nNUC*2**TAUY 

CALL PM (EMI ( \ ) # DNuBfEMx*Gl ( 1 ) ) 

CALL fM3d<EMX *£MR»XS*G1 U) > 

PBsP? (EMR) 

OUH- { p A * PR ) * ( Zr^Za* /4 * 
r>XLFs»DUM#(XB-XA ) 
nXTHsnUM#ABS { yr-ya) 

YRB= (yA+YB) /2*~YSHF 

XBR = <*» ( (XA + XS)/?a**XSHF) 

dxmom-^yrb^dxtm+xrb^dxlf 

XLIF=yLTF+DXLF 

XTH“XtH+dXTH 

XMPMsxMOM*DXMOm 

XAsXB i Y A = Y 8 % 7A-ZB $ PA^PB 

FMAsFmB 

PX*Pa*PIWF 
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PX*PX/PJNF 

IF (IPp.FO* 1 > GO TO 5b 

WRITE (6* 100) XA* VA»ZA # Px*EMA,XLlFt XTH,XMOM 
55 CONTINUE 
40 CONTimUE 
60 CONTINUE 

xth=xth+xthv 

XLTF=xLIF*XLV 
XMOMsyMOM+XMOMV 
WRTTF (6 *4203 ) xLTF 

4203 FORMat</// 20X*TOTAL LIFT a*El2,4) 
wRITF( 6, 1 0 08) y TM 

100B FORMAT </20**TH&UST (FROM COWL AND VEHICLE UNDERSURF AC?) a #El2t4> 
write ( ft t 4206) XMOM 

4205 FORMAT</5?0X*MOmENT (FROM COWL AND VEHICLE UNDERSURFACE) =#El2 8 4) 

IF ( T H x • L T • 1 • E ♦ r» 9 ) GO TO 4201 
WRITE (6*4200) 

4200 FORMAT (/20X*EXIT PROFILE AND TmRUST AnD MOMENT (FROM SIDE wALLS) N 

iot calculated#*) 

C,0 TO 42)2 

4201 IF (THx,GT.O. ) f?0 TO 4205 
WRITE (6*4204) THX 

4204 format (/20X#ThrUST (FROM SIDE WALLS) =*E 12 #4) 

PMOMi „0 # 

WR J TE ( 6 ? 42 1 1 ) PMOM1 

4211 FORMAT (/20X*MQM£NT (FROM SIqe WALLS) ^E 12 #4> 

GO TO 4212 

42qB 0 THGstHX-XTH 

WRITE (6*4204) dTHO 
XTH*ThX 

PMOMi*PMOM-XMOM 
WRITE ( 6 f 4211) PMOM1 

4212 CONTIMUE 

PTHL® (THRMAX-XTH) /THRMAX*1()0* 

WRJTf(60563) pThL 

1563 FORMAT (/40**PERCENT THRUST LOSS s*F6,3) 
return 
END 
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SUBROUTINE THM ( E M ! ? P 1 * G ? P I NF « Y V 1 9 YC 1 1 D Y V ? X V 7 9 T H RM A X ? X c 2 ) 
COMM0N/D3/AZ f B7*CZ 
?aAZ»xV2<H*2*B2*XV2 + CZ 

a1«*yvi-yc1 

Z2=Az*XC2*XC2+PZ*XC2*CZ 
A2 3 DYv# <7+Z2 > /?« 

FMl=l .+ (G-l • ) /2**FMl*#2 

OUM=<;>*#fM1/<G*1 # > )**( (G + X,)/2 f /«3-U) ) 

A ST 3 A ]. *E'M 1 /DUM 
PTOT = pl/FMl**(r5/(l *-G) ) 

EM2*EMi#S^«T ( AP/A1 ) 

IF (EMi.GT.3 * ) Em2«(A2/A1 ) 

AF K A ? /A ST 
X TM = 0 
• 10 CONTINUE 

FM?=1 .«■ {G-l* } /?c4*EM2**2 

OUM= ( p**FM2/ (6 + 1 * ) ) ## ( <G + 1 « ) /2, / (G-l • ) ) 

aft=om m ^fmh 

ERAs(aFt_AF)/AF 

IF(AB5(ERA) .LT,l.E-03) GO TO 2 n 

CALL pRROR(2000tITMfEM2»ERA*l.l*EM2l»ERAl) 

GO TO 10 
20 CONTINUE 

FisP U*G^EMx #tf 2) 

PHsFM?*# (g/ ( I .-G)) *PTOT 
F2=P?*A2«- ( 1 „ +G^FM2##2) 

THRMAx 3 F2-F‘1«PINF^ (A2-A1 } 

RpTURM 

FND 
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SUBROUTINE XTH(XP,Y2*TH2,DNUF* FMl 9 G1 , pTQT * TTOT ♦ THS2 * 
iTHXl , XMASS1 *THX) 

n I mens ION XNU (50 ) fX(50)vY(50)«?(50)*EM(50)9P(50)fTH(56)9T<50) 
COMMON/P /PINE 
C0MM0 n/D 3/ a Z »87»CZ 
COMMOn/5HF/XShf* YSHF 
COMMOm/PmPM/ PmOm 

nT mens I on DNUE(IO) »X2 (10) * Y2 { in) * TH2 ( 1 0 ) f ( 03) »Gl ( 03 ) 

GETZ (y5) s AZ*X5*X5*0Z*x5*CZ 

THX-n * 

XF(AZ^HZ,EQ.O,> return 

INDEyaO 
180 CONTINUE 

DO 1 0 I a 1 9 5 

J 2 s I w S 

IF(I # £G.5> J2 = ?6 

PiO 2ft J = jlfU2 

IF ( InoEx.EQ,!) GO TO 86 

RAT*F| t OaT < J-Jl ) / FLOAT ( J2«Jl ) 

XNU ( J)«DNUE ( l ) oRaT* <DNUE ( 1*1} ~n>NU£ U > ) 

X<J)s X 2 n> *RAT»<x2<I*l)«>x2(l> ) 

Y ( J) =Y2 ( T ) ♦PAT* (YR ( 1*1 ) *Y2 ( I) ) 

TH ( J)*TH? ( I ) ♦RATO (TH2 ( I ♦ 1 ) «TH2 tl)) 

CALL pM(EMl (1) «XNU<J> t EM ( J ) , ($ 1 f 1 ) ) 

CALL fM3o(EM( J) 9fMX»X(J) »G1 (1)') 

FM< J) = EMy 
7 ( J) sc.ET? (X ( J) ) 

8G CONTINUE 

DUM-] ,*(G1 <l)-l . )/2.*EM<j)#EM(j) 

DG = G? f 1 ) / < G 1 (!)■•],) 

P ( J) -pTOT/DUMttttOG 
T (J)jsTTOT/DUM 
20 rONTTNUr 

10 continue 

HUMsGj ( j )**2*SIN(THS2 ^Tw (j) )/S lN(THS?) 

FlsP(l ) *pM ( 1 ) /SORT (T ( 1 ) ) *SlN (TwSZ-TH ( 1 ) )/ 

1SIN(TH52)*Z (1) 

U"1 * wplNF/P ( l ) 

F 1 “P (1)*(U ♦oUH^cOS ( TH <1>>)#7 (l) 

Wl=P ( 1 ) * ( U/TAN { THS2) *DUM^SIN ( TH < 1 > 

XMAS$2“0 * 

THX2 = ft o 
PMOM = ft # 

DO 15n J*2*26 

p2»P ( I ) *FM ( I ) /«;QRT (T ( I ) > * 

1SIN(ThS£_TM (I) )/srN(THS2)«z (T) 
nUM=0l ( 1) ( T ) #*2 *SIN (THSs-Tm < I > > /SIN (Th$2) 

Usl,-plNF/P U) 

F2 = P (I)*(U *0ijM*C0S(TH (I)))*7 (I) 

W?=P(t> » (U/lAN (THS?) ♦OUM*Sl N ( T m ( I ) ) )*7<I) 

THX2 = tHXP« (FI + p2) * ( Y { I ) »Y { j ■» 1 ) ) 

XMASS2=X M AS52+ B 5ii‘{El^F2)*(Y ( I ) Y ( I ** 1 } } 

PUMls.5* (F1*<Y { I~1 >-YSHF) *F2 <My ( X) - YShF) ) 

OUM2b. 5* (wl*(X( I-l ) -XSHF) *W2 *(y (11-XSmF) > 

PMOMkPMON^ (DUM i -0UM2 ) *(Y(I)-y<!-1>> 
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E1=E? $ F1=E2 % wi=W2 
150 CONTINUE 
100 FORMAT (7E13*5) 

rATMsxMAsSH/XMASSI 

EMSsRaTm«1 

IF(ABS<EMS).LT,l.Er03)60 TO 46?5 
TNnEX»lNr)EX*l 
IF(INoEx.GT. 1)60 TO 4625 
DO 170 1=1*26 
FMl=FM<FM (I)fGl(l)) 

FM2=FMl/RATM 

TTM=3 

FMT = RaTm<*EM (I) 

19 7 CONTINUE 

FMTssFm (EMT*G1 ( 1 ) ) 

PPM=(fM?-FMT>/fM? 

IF ( A8S (£RM> *LT* 1 oE~03) GO TO 1 7 f 
DUMD-i ,/RATM 

CALL pRROR <6f ITM* EMT* ERM*DUMd*eMTi »ERM1 ) 
GO TO 197 
171 CONTINUE 
EM (n*EMT 
170 cONTTmUe 
GO TO 180 
46?5 CONTINUE 

THX = ThJ< 2'-THX i 

RETURN 

FND 
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SUBROUTINE £m3d <EM,£MX*X*6) 

C0MM0N/D3/ A^*r/«>c2 
Z = A 2^X*X+t?/:*X + c Z 
fiN= ( 54 1 • ) / ( G* i • ) / 2 * 

Ys() ,* (G-l • ) **GN 
Y=Y/ E M/ { (g*1 • ) /2. ) ^*gM 
A«Y*Z 

fr M X = £ M ** S Q R 1 ( Z ) 

I T M = 1 

10 AT« ( 3 .♦ (G-l • ) /P.*EMX*eMX) ##GN 
ATsAT/fcMX/ ( <g + 1 • ) /2- ) *#QN 
ERMs ( 6 T- A ) /A 

IF (ARs’tERM) • L r * 1 .E-03) 60 TO 2 ^ 

CALL FRR()RUu» Itm» £mX tE«M* 1 * 1 ?EmX1 jErmI j 
00 TO lo 
20 CONTINUE 
RETURN 
FND 
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SUBROUTINE fix f Xl ,V 1 ,THi ,XMU 1 » ¥2 * YZ f THZ • XMU? »X 3 t Y 3 * TH3 , XMU 3 »EM 3 f G > 
SL 1 = T A N { TH^Xmi.! 1 > 

Sl2 ss T a N{ThZ-Xm»»2) 

CALL GfcH ( XI f Y l , SL1 * X2 f Y2 # SL2 , X 3 , Y3 ) 

I T M s 0 

10 rONTlhlUr 

CALL £ m 3 0 < t M 3 9 E M 3 N 5 X 3 > G ) 

X M U 3 = a S 1 i\i ( 1 • / £ m 3 M ) 

SL 1 3 = SL I /2 • ♦ • b tfT A N < TH3 + XMU3 ) 

jP(Xmu1*E'G*0*)S-L13«Sl1 

5L23-$L2/2 . + . AM ( [ H3«XM^3 ) 

r p ( X t j 2 * £ Q • 0 • ) s £ 23”S[ 2 

CALL gLM ( x 1 9 Yi ,SLI3 ,x2,Y2,SL23 s x 3, y3) 

jp ( j T m*E&* 1 ) T 0 ^0 

T T M “ 1 

GO To 10 

20 FH 3 = EM 3 M 
Pp TURN 
END 
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SUBROUTINE Rtf { EMI * DNIJ t Ert2 *G > 

GG=SGRT< (G*1.)/(G-1.) > 

Xtf lr s O K T (Etf * ) 

XNU 1 sGG#Al'AN ( AM} /GG ) -A TAN ( XM\ ) 
fMP^qNjU/ { J. • 5-XmU 1 ) ~em1 ) +emi 
X T 3 - 1 

10 Xtf2=SfJ«T * ) 

[ )NUT = rG*(A t aN(VM2/gG^-aTaN(XM1/ 6 g) > + ATaN(XM] ) ™aT,aN ( Xtf 2 > 
E k N U - u NU „ 0 N u T 

IF < A BS <F«NU) ,tT • 1 *£-04) qQ Tq 2ft 

CALL ERROK (3, U3,EM* f E«NU, 1 , 1 1 ,EM2l *EPNU1 ) 

60 To 1ft 

20 CONTINUE 
RETURN 
END 
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SUBROUTINE POLY (XVI * YV 1 * TH I t X. V 3 ? Th V3 » XV2»THV2» 
lYclfYcl^C^’Trir 2 * x lJ 

COMMOM/y/Al 1 j A ] ,2 9 A i 3 t A 2 1 9 A22 * A?3 * A3 1 » A J£ * A3 3 
cOhmOn/POl/ iPr, rPOuy 
n I ME Ni 5 XO Ki * U 0 ) * Y ( 1 0 > , yH ( 1 0 ) 

IF ( iPoLY.gO. 1) qO TO 20 
I T s 0 

X (1> = X v i 

X ( 2 ) RXV1+*W1* (XVP-XV1) 


X (3) =xY3 
X ( 4) -X v 2 


T H ( 1 ) = T H I 
TH.( 2 ) sTHV 3 


Th (3) =t h v3 
TH (4) bTHV2 


Th ( 1 ) =-th t 1 ) 

Y ( 1 ) ssYCl 
Ns 3 


$ Th(2)*-Th<2) S Th(3)=-TH(3) 


$ TH( 4 )=-T H (^) 


WrItE f b* 1 00) 

100 FORMAT (///30X-VEHICLE COOHO I NATES*/ 

1 1 5x*X tt ^X # T^ #4 X* A* 15 x*COOrD INAtEs** ) 

30 DO JO Isl*ISi 
X? = X ( 1*1 ) "X 1 
X I = X ( j ) - X I 

Aa ( T Ah! ( TH ( l + 1 ) ) -Tan ( Th ( I ) ) ) /? f / ( x2-X 1 ) 
RbTAN (TH ( 1 + 1 ) ) -2,*A*X2 

Y ( *+ 1 > =A#x2*x2 + S* x 2+C 
WRITE (&* SW) A(I)*KU + 1)»A»B*C 
6* p 0 R M A T ( 1 QXb^l 1,3) 

10 continue 

IF(IT,Eq m )GO tu 20 
IT-1 

X (1 ) rxCl^X(2)=x U) +*01* (XV2-XV1 ) 
y (3) =xc2 
TH< 2) = V H C 2 
Y H ( 3 ) a j H ( 2 ) 

t H(2)=- t H(2) $ t H(3)s-jH{3) 

Y(1)=YV1 

Ms? 

wRJTe ( b, 101 ) 

1 o’ i FORMAT (/3 0X#COV'L COURojNaTeS*/ 

115X*X#4XttTU#4X*X*l5X»C00K0INATFS*) 
fiO TO 3r 
2U CONTINUE 
RETURN 
F N 0 
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SUBROUTINE GEM (XA * YA * SLA* XB* YB* SLH »XC, YC) 
xCs { Yf3-YA^i>LA*XA-SL Utt Xa) / < Si A-SlB) 
YC=YA*SLA* (XC-xA) 

RETURN 

END 
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SUBROUTINE ERR OR (ItIT*X»ER*F*v. j. igRl) 

J T = I T 4. 1 

IFflT.LT.15) GO TO 
WRITE (6 m 3 ) 

1 3 FORMAT (*fRKOR T£ST NUMBER 
WRITE (6,20) I 
20 FORMAT (15) 

STOP 

1 2 IF(lT*GT t2 ) GO TO j4 
g;Rl=£R 
X lsX 

X *= X * F 

T F { X . F U * X 1 > X=X + .l)2 

RETURN 

14 XDsXl-ERl* <X-X1 J /<ER-ER1) 

Er1*Er 

XI sX 
X«X D 

return 

end 
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PROGRAM nOZD ( INPUT 9 OUTPUT »TAPEr= INPUT, TAPE fesOUTPUT) 

COMMON/ A/ PItTNltEMlfGl*Tl*HX»OHl*UI 
COMMOn/B/XA,YA, XR, YH,XC 9YC« Xn* YD» 

ITHA »ThB9THC9TH ! S9EMA»EMb»EMC9EM^ 9 XMUA*xMU8»XMUc»XMUD 

COMMOM/C/DNUV9DNuC?XC 1 * YCl 9 XV 1 1 YV 1 » ONi jL 
COMMOM/D/ PHI 

COMMON/BP/^A9Pn»PCf p D*TAf Tg • TC , TD o G a t r,8 ♦ GC t GD t UA »UQ , Uc t UO #«H A * 

1 RHR 9 RMC » PHD 
COMMOnj/D?/ &ZfPZ*CZ 
common/v/aV, av,cv 
commonj/p/pinf 

COMMON/SHF/XSHFf VSHF 

commom/pol/ipr, I poly 

COMMOn/ThRMaX/THRMAX 

DIMENSION P2 ( 1 0 > a TH2 (10) ,EM?(ln> , TZ ( 10 ) t GZ ( lO ) 9X2 ( 10 > , Y2 ( 10) 

1 , 2 ? (in) , R H 2 ( 10) 9 U2 ( 1 0 ) 
n I MENS I ON HOL (11) , HOLE (6) » DNUE (10) 

DATA HOl E/2HE0 f 2HEl t2HE2,2HE3 9?HE4*2HEb/ 

DATA M 0L/2HVl92HClf2HA *2HR ,2 hC ,2HP t 2H V3 , 2 hC3 • 2hF ,2HG ,2HV4/ 

GFTz ( X 5) -Az^xS^xR+B2 # x5+Cz 
96 FORMAT (7El3*5) 

99 FORMaTl/i 
100 FORMAT <flFlO,0) 

READ (s, 100) PI ,THI ,EMI,II , PH! , PINE 
H I sF H ( P I ,PHI ,TT) 

GIbFGaM (T l,PI,PHI) 

RHI=RHEO (HI ,PI ,PHI » TI ) 

ASSORT (gI^PI/RHI) 

UI«EMi*aI 

TaURc«°* 

RE AD ( r j 1 O0 ) XVl’YVl’XV^fX^l^YCl^XCSfDYV^UNtjC 
0NUC=dNUC/57.3 

READ(r* 5r 21) IcFtOXC* JTF*DTH*IPRfIPOLY 

5921 FORMAT (I59ElO,09t5tElO a O,2I5) 

PFAO(^»loO> XSHFtYSHF 

READ (5* 100) AZ,8Z,c 2 
ZClaGETZ(XCl) 

7V1 »GF T Z (XVI ) 

7 C j? = Gf*TZ(XC2) 

WrITE(6,b922) 

5922 FORMAT ( 1 H 1 / ////) 

WRITE (6, 1020) 

1020 FORMAT ( 14x#E QUILIBRIUM FLOW NOZZLE D E 

1 S I G N <f ) 

WRITF(6,400) Pl,THl cEMI »TI,PHI ,PINF 
400 FORMAT ( ///30X»*INITIAL PROr I LE <V/3 X*PPEsSURE *6 x#ThETA<*9 x«MACH# 
14X*TEmPERATURE#6X*PHI»9X*PTNF*/6Ei3#5) 

WRITE (6,5930) XVI , YV1 »XV2*Xcl »YCl *XC2,OYV 
593 0 FORMAT (///f>X#xVl*l 0X*YV1 #1 0x*XV2*l OX*xCl*lOx*YCl*lOX*xC2*I ox *O yV*/ 
17E13.P///) 

WrIt^ (6*5923) xS WF,, YsHp 

5923 Format (^X^mOmEnT AX.ts*/12X**x =*E13*5/i2X*Y s#E13„5///) 

WRITE (6,1001) AZ,0Z,CZ 

1001 fORM a t (qX^LaTePaL eXP A NSlON e q HaTION^/9X»2?mZ(X) = A Z#X»*2+0ZftX+ c 2 
1/12X*AZ s # EI3,5/12x # Bz s*E13 # 5/12x*Cz s*E13 5) 

‘ XC 22 skC 2 



TR 188 


Page A3-2 


THCl=rHl "*DNuC 
YC2»Yc^TAN(THC1 )* JXC2-XC1) 

YVP=YC2<frDYV 

ONUVsaTAM ( (YVg-YVl) / (XV 2-XV1 ) ) «THI 

no SOOO T xc- 1 « I CF 

XCP = yC2? + float (IxC-^1 >*dxc 

CALL THM (RHI ,UT »PI *PINF f YVI *YC| *DYV* XV2 t THRMAX *ONUV tONUC f XC2 > 
WRITE (6*6364) tHPMAX 

636 A FORMATUHl31x*T0EAL THRUST =*£12,4) 

YC?=Yrl*TAN(THCl ) <MXC?*XC1 ) 

YVP“Y C 2+ D YV 

no e>ooo jtsi,jtF 

THX=] .E+10 

THJTsFLO'AT (JT-.1 ) *OTH 

nNUV=A T AM ( (YV2-.YV1 > / ( XV2-XV1 ) ) -THI *TH jT/57* 3 
CALL pM { nNUC*PCl * TCI sFMCl ®GC1 ?!!C1 ?RHCl ) 

CALL pM(ONUV»PV 1 *TV1 tEMVl fGVltllVi*RHVi ) 

TT?=1 

1 0 0 0 continue 

CALL cHAP 
DNUArO, 

ONUL)ar>NUV*DNUC 

7A=GET2{XA> 

7fi=GETZ ( XB) 

7c=fipT2 (XC> 

ZD=GetZ(XO) 

pat=o, 

THC3RT HC1 

CALL pM(DNUCtPC3tTC39EMC3tGC3»MC3*RHC3) 

CALL FIX (XC1 »Ycl *THC1 * O, »X« • YB *TH8 tXMUfi * XC3 » YC3 * THc3 9 XMUC3 e £MC3 « 
1GC3* 

1PC39TC3* PHC3*UC3,DNUC) 

7C3*GFTZ<XC3) 

T T 1 7 6 0 

jF (XC3 *l t » XC2 > GO TO 177 
IT176 s 1 
TAUCeOt 
X.C3 = xC2 
ONUEObONUC 
GO To 1 78 
177 CONTINUE 
IT9»1 
RA T“ o 1 

176 CONTINUE 

XMUBf>sXMUB*RAT» (XmUD"XMUb> 

THRD = tHPURAT* (THP“THB) 

XT~Xp ♦ P AT* < XQ=»XB ) 

YT=YR*RAT* ( YD-YB) 

TAUC=rAT*ONUV 

ONUWsnNUC*T A UC 

CALL pM (ONUW f Pw,TW«EMW*GW*UW f RwW) 

THW = ThC3-*-TAUC 

CALL FIX (XC3 f YC3»THC3,XMUC3f XT*YT t THBOtXMUBn«XW*YW»THw»XMUWfEMW 9 

l.GW« 

lPW^TW.RHWtUWfOMUW) 

' DNUE0aOfMuC+2# *TAUC 
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178 CONTINUE 
THC2=tHCi 
THEOstHc" 

CALL nM ( DNUEO»PEO,TEO,EMEO,GEO f uEO»RHf;0) 

7E0«GETZ(XC2) 

RQ3D = pHE0<*UE0/7E0 
DNU3e s DNUE0*SQrT(2E0) 

ITEO=i 

200 CONTINUE 

CALL pM(DNU3E*PE0 9 TEQ,£Me;0,GFO,uEO*RHfO> 

P03T=pH£n^uEO 
EROEs ( RQ3T-R03H) /RQ30 
IP (AGS (£RQE) • LT*1.E‘"03) GO TO ?l0 
CALL ERROR (200, I TEO * DNU3E , ERQE , 1, 1 tDNU3El .ERQEl ) 

GO To 200 
210 continue 

XMUEO-ASIN ( 1 ,/EMEO) 

IF(IT|76.EQ,x) GO TO i75 

SL w C2s; (TaN(ThW-XmuW) *TaN ( The0~xmU£0 ) ) /Z • 

VC2T-YW+SLWC2* (XC2-*XW) 

FR9=YC2-YC2T 

IF(A8 s (Er9).Lt,1.E-03) Gq TO 1?5 

CALL F R R0R<4*lT9»RAT*ER9a,n*PATltERRl> 

GO TO 176 
178 CONTINUE 

THV3 = tHI +ON1.IV 

CALL pM (f)NUVtPV3*TV3*ENV3tGV3tiiV3»RHV3) 

THV 1 “THV3 

CALL FlX(XC ? YC,THC» XMUC, xVX 9 yV 1 ,TMV1 »0. tXV3*YV3,THV3»xMUV3,EMV3 # 

1 G V 3 1 

1PV39 Tv3*RHV3i U v39f)NUV} 

AV-YV3 % BV-TAN(THV3) 

CV= <YV2-aV- 8V<* (X.Vg-XV3) ) / (XVg-XVa) **£ 

7V3=CrTZ (XV3) 

THV2aATAM(BV^2 # #cV*(XV2-xV3) ) 

IF ( TH\/2oLT * o • ) GO TO 462 
TAURV»Thv2-THV3 

T T 7 ® 1 
I T A U ts h 
rat=i , 

54b CONTINUE 
530 CONTINUE 

XZ=Xc+RAT* (XQ-xC ) 
yz»yc*Rat + m YQ^yC) 

TAU«PAT*r>NUC 
THZ=ThC+RAT* (TnD^THC) 

EM7 b EMC*PAT* (£MD~E m C ) 

XMUZ = /\SlM ( 1 «/£MZ ) 
t Hy=t h V3-t a U 

pNUY= n NUV+TAU 

CALL pM{DNUVfPYtTY*EMYfGY9UY,RMY) 

CALL f IX (XZ*YZ f TwZ*XMUZ*XV3»YV3*T H V3*XMUV3*XY*YYtTHY»XMUYt ? mY*gV* 

1 P y * T y * RH y 9 U y ♦ DNU y ) 

ONUE5ssONUV*2*#tAU4TAURV 

Tnf5=THV? 

• CALL PM ( DNUE5 t PE5*TE5fEME5,GE5,UE5»RHFS) 
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7E5*GETZ(XV2) 

R03d=rHE 5#UE5/7£5 
nNU35aDNUES*SQRT ( Z E5) 

ITES=i 

230 CONTINUE 

CALL pM(0NU35tPE5t TE5*EME5*GE5 f uE5»RHF5> 

RG3TspHE5#UE5 

rR05= (RQ3T-RQ3n) /RQ3 d 

TF(ABS(EPO5)*LT*l.E-»03) GO TO ?30 

CALL r R ^OR (220o ITE5*DNU3BfER05?i« 1 *0NU35l »EP051 ) 

Go To 230 
.330 CONTINUE 

XMUE5=A$T.N <1 •/FMEB) 

SL2Tsb (Tan (THY*XMUY) *TAN(THf5*Xmue 5) >/2» 
yV2T=yY+SL2T« (XV?«XY) 

IF(IT7.EQ*1*AN0*YV2T«GT*YV2) Go TO 54? 
pZ=Yv?T-YV2 

IF<ARS(F7) *LT t j # F„03>GO TO 540 

IF (ONuC „ EQ • G • ) GO TO 540 

CALL ERROR (5# ItZ#RAT»FZ» , 1 #R«Tl*Ezl) 

00 TO 530 
542 ITAU = i 
540 CONTINUE 

DNUF psDNUV *DNUc 

CALL ■pM(DNUE2fPE2*TE2«EME2.GE2*UE2»RNF2) 

THE2stHI ♦QNUV-oNUc 
THFsThE? 

RHF»RHE2 % UF=UE? 

CALL FIX (XC3» YC 3 ♦ THC3 • XMUC 3 * XD i YD * THD * XMUD * xF * YF ♦ THF * XMUF 9 £MF 9 
IGF • 

lPF*TF t RHF»UF*DNUE2) 

SLE= (YV2-YC2) /(XV3-XC3) 

THE s AT‘AN ( sle ) 

UFT=Uf 2 * RHFTsRhe2 

CALL plx (XC3tYC3ffHC3 t xMUC3txC2»YC2tTHE,0 # .xFTtYFT*THF f xMUFTfEMFT, 
1GFT, 

^ppT^ T fT * RHFT f U fT »oNUE2^ 

TF(Xf.GT.XFT) c,0 to 300 

CALL Fix (XF,YF.THF * XMUF tX C2*yC2 f THE, 0, »XE2,YE2*THE2 f XMUE2,EME2, 

1 G E 2 • 

lPE2» TF 2 ? RHE29UF2 t f)NUE2) 

7F=Gf t 2 ( xF) 

SIF s TaN(thF-XMhF) 

CALL fiEM(XFtYFt5LF*Xc2tY C 2f SLE*XE1*YE1) 

ITElen 

36B cONTimUe 

RAT=(YEl-YC2)/(YE3-YC2) 

ONUEIbDNUV+DNUC* ( U-RAT)*TAIJC 

CALL pM { DNUEl ,PEl,TEl ♦EMEUGEl ,UEl#RHrl ) 

THE1=tHI ♦ONI.IV-ONUC- ( 1 ,-RAT) #TAUC 

CALL FIX (XC2 • YC2» THE • 0* *XF*YF*THF * XMUF »XE1 tYEl • THE 1 fXnUEl •EMEl • 

IGF 1 « 

1 PE 1 ?TFl*PHElfUEl oONUEl ) 

TF ( I TFl *FQ. 1 ) GO TO 369 

ITElrl 

GO TO 368 
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369 continue: 

CO To 3m 

300 CONTINUE 

CALL FIX (XC2i YC2*THEt oo *XD»VD»THDtXMUo*XE2*YEa#THE2»XMUE2?eNE2t 
1GE2* 

lpE?*TE2tPHE2*uE2*ONUE2) 

SLCL' = TAN (THC3 + YMUC3) 

CALL GEM (XC3»YC3tSLCE#XC2tYC2*«;LEfXEl»YEl> 

TTflsn 

301 CONTINUE 

RAT = ( YE1 -YC2 ) / ( YFP-YC2 > 
nNUEl=ONUC*TAUC^PAT*(DNUVTAUC) 

CALL pM(DNUE 1, PE1*TE1,EME1 *GEl,uEl*RHcl) 

THE 1 . =TH I •*ONUC*TAUC^RAT*{DNUv^TAUC) 

CALL FlX(XC3tYc3«THC3tXMuC3*XC?»Yc2»THE»0, f xEl»YEl*THpl*XMUEl*EMEl 
1 * GE1 * 

1PE1 tTEl»RHEl»UEl«>DNUEi ) 

TFUTFUFQtl) GO TO 300 

TTElsi 

GO TO 301. 

380 CONTINUE 
3 j 0 CONTINUE 

THE3=tHf? 

DNUE^ssONUEZ 

CaLL pM (f)NUE3f PE3 9TE3 *EMe3»GE3^UE3»RHe3) 

XMUfSsASTN ( 1 */fMf3) ' 

IF(ITAU # NE,1) GO TO 311 

X G = X Y $ YG h YY % THG-THY $ RHG^RHY $ Uo s UY $ £MG=EMY $ XMUG=XMUY 

pg=py $ tg“Ty 
xgt=Xc, 
go To 35*8 
311 CONTINUE 
THG=THE3 

RHG=RhE3 S UG-UE3 

CALL pIX <Xo»Yo*THD*XMUD*XV3tYV3,THV3»xMUV3tXOtYG»THG*XMUGVFRG« 

1GG, 

1 PG f TG ♦ RHG * UG * DNUE3 ) 
rHGTsrHe3 $ IJGT=UE 3 

CALL FIX (XC2*YC2tTHE9 0. » xV3*YV3*THV3,yMUV3, j(GTtYGT*THG»XMUGT9EMGT9 

1GGT* 

1 PGT • TgT* RHGT»UgT$dNUE3) 

IF <XG # GT # xGT ) GO TO 350 
3589 CONTINUE 

CALL FIX (XC2* YC29THE*o* •XG*YG«THG*XMUG*XE 3»vE3tTHE3*XMUE3*FME3* 

I G F 3 9 

1PE39TE3* RHE3»UE3 f 0NUE3) 

ZGsGFTZUG) 

IF(YE3«LE*YE2) GO TO 35g 
SLEA-tAn (THG+XMUG) 

CALL fiEM (XG* YGfSLE4*XC2*YC?tSLF*XE4tYEA) 

T T A~ 0 

350 CONTINUE 

RAT= (yL'a-W?) / ( YF3-.YV?) 

r)NUE AsDNtjV > qNUC f 1 .-RaT) * (T A URv*TaU) 

CALL PM (0NUE4i P£4 * TE4 *EME4* GF4 ,uE4t RHF4) 

THF4s:tHI ♦ ONlJV-DNUC* ( 1 . -P A T ) * ( T AUR V ♦ T AU ) 
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CALL FIX (XG* YG* THG* XMUGf xC2f YCp.THE * n . t XE4 ♦ YE4 * THE4 » XMUE4 » EME4 , 
lC£4« 

1P£4* 7F4,PHE4,UE4 f HNUE4> 

TF ( IT4*E0. 1 ) GO TO 359 
T T 4 ” 1 
GO TO 358 

359 CONTINUE 
IF(lTflU f r0.o) GO TO 360 
SIV3?- (YV2-YV3) / (XV2-XV3) 

SLG»TaN(THG + xmijG> 

CALL rEm ( X V 3 1 YV 3 • SL V3g * Xg * YG • StG • X V4 • YV 4 ) 

TTV4 = i 

362 yV4eAV*PV<* (XV4 W XV3) *CV* (XV4_XV?) <>*2 
THV4 = /\ T AN <8V*2,^CV» <XV4**XV3) ) 

TAUPV4=THV4-THv3 

PNUV4bDNUV*2.*TaU*TaURV4 

CALL pM (DNUV4, PV4 t TV4*EMv4,GV4aJV4»RHv4) 

CALL FIX {XV4 »YV4? « 7854*0. »XV4»YV4tO* »XV4*YV4fTHV4t xMUV4*FMV4» 
1GV4,PV4, TV4, RHV4,UV4 9 DNUV4) 

YV4T”y^+«5' < *(TAM( tNG*xNUG) + TAN { TH\/4*XMnV4 ) } *<XV 4 ~XG> 

FRV4= ( YV4T«*YV4) /YV? 

IF (ARS (£RV 4 ) *LT*1 *E~03) GO TO 36 1 

CALL FRFOR(36l,ITV4,XV4 f ERV4,l.l9XV41,ERV41 , ) 

7V4=gfTt: { x V 4 ) 

IF ( I T\/4 e pQ . I ) R 0 TO 361 
GO TO 362 
361 0NUE4s0NUV4 

IF (YE 3 *LFe Y e2) GO TO 460 

CALL pH(0NUE4,PE4 f TE4tEMe4*GE4,uE4.RHc4) 

THE4=tHV4 

CALL FIX ( XC2 * Yq 2 * THE * 0 • *XV4 t YV4 ♦ THV4» xMUV 4« XE^ » Y^4 * T Hfr4 # XMUE 4 9 
1 FME4 , G-4 ? 

lPE4»TF4 9 RHE4tUE4«DNUE4) 

GO To 360 
350 CONTINUE 
tT4 = 0 

CALL FIX f XD t YP,THD f XMUD f XC2 t YC2 t THE 9 0. # XE3 f YE3,THE3,XMUE3 f EME3 9 
1GE3* 

1PE3* Tf3 ? RHE3*UE3»ONUE3 > 

JF ( Yf3»LF« YE2) GO TO 460 
5LE4 a TAN <THV3~XMUV3) 

CALL oEM <XV 3 * YV 3 *SLE4 * XC 2 ♦ YC? * C LE 9 XE 4 * Y£4 > 

371 cONTimUe 

P^t= ( V E4«yV2) /(YE3-YV2) 

DNUEa-DNUV+TAU+RAT* (DNUC-TftU) 

CALL pN (oNUf.4 n PE 4 ? TE4 . pMfr4 9 gE4 o(JE4 * RHf4 ) 

THE4 c tHi' ♦ONljV^TAU-RATi < DNUC^TAU) 

CALL FIX (XC 2 tYC 2 »THE *0* * X V 3 * y V 3 * TH V 3 * XM U V 3 * x E 4 * Y £ 4 * THp 4 * X MUE 4 1 E ME 4 

1 AGE 4 * 

1»F4 * Tr4 p RHEA 9 UE4 ,0NUE4 ) 

IF (TT 4 .EO. 1 )GO TO 3?0 
T T 4 “ ? 

00 Tp 371 
370 CONTtn|UF 

360 CONT’mUf 

TF (YFlLE.YE?) GO TO 460 
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X.?a>**C? s X2(2)sX£l $ X2<3)=VF2 % XM4)=Xr3 * X2(5>^XE^ 

v?(6; 2 X>'? 

vp(l) r YCP S Y 2 (?) = V E 1. % Y2(3)=:vP? $ Y?(4):sYF3 $ Y2(5 )-yE4 
Y2(6) s YV? $ EM^(i)=EMEO $ EM ? ( ? ) =fME i $ E^2 ( 3 ) “EMEE S £M2U) 3 E w E3 

£M2(«;)*EME* $ FM2(6)«EME5 $ TH?(1)=THfO $ TH2(2)sTHrl 

TH2(3)=THE2 $ TH?<4)=:THE3 $ TH?<5)=THf4 $ TH2(6)=THE5 

DNUE ( ) ) =qNUEO * DNUE< 2> *DNuei * DNUE ( 3 ) -DNUF2 S DNUE ( 4 ) =DNUE3 
ONUE (5) rDNUE* $ DMUE (6) spNUEB 

PH2(1)=RHE0 $ RH2(2)=PHEl S PH?(3 )=RHf2 $ r h 2 ( 4 ) =;RhF 3 
RH2 { 5 ) sRHE4 $ rH?(6)s:RHE5 

U2(l}=Ueo S U2 (2 ) “UE 1 $ U2 ( 3 ) -I.JF2 $ U?(4>»UF3 $ U2<R)*UE4 
U2(6) s UE5 

G2{1)s;Geo $ G2 <2) esGEl $ G2{3)=CF2 $ G?(4)s:GF3 $ G2<5)=GE4 
G ?(6 )s3E5 

P2U)=P£0 S P2f2)=PEl $ P2(3)=»E2 5 P?(4)=PF3 $ P2(5)=PE4 
P2 (6) =PEB 

T2(1)=TE0 $ T2 { 2 ) = TE 1 $ T2(3)=TE2 $ T?<4)=Tf3 $ T2(5)=TE4 
T 2 ( 6 )s TEB 
TMDEX.O 
18° CONTINUE 

DO I60 I si *6 

72(I)=A2^X2(I)ttH>?*PZ^X2(I}^CZ 
160 C ONTinjUe 

THS1 = i ,5707963 

TF(XVi#NE*XCl> THS l=ATAN ( ( Y V \ -YC 1 ) / ( Xv 1 -XC 1 ) ) 

THS2= A TaN{ (YV2-Yc2)/(XV2-X C 2) ) 

7lGs,e;*(GETz <xCl> *GET7 (xv* > > 

XMaSSi =RHI # UI*SIN (THSi-ThI ) /S I N< ThS 1 ) # Z I Q» ( YV 1 - Yc X ) 

THXl= (Rhi # UI«ui*sin<ThS1-ThI) V os ( THI ) * (PI-PINF> *SlN <THSl ) ) # Z I G* { Y 
lyl-YCl )/SlN(TH<;l) 

CALL XTH(X2»Y2, TH 2 fONUEtPlNF t TwS 2 tTHXi ? XMASS 1 9 THX ) 

TSH=0 
GO TO 461 

460 CONTINUE 
lSH = j 

wRITf ( 6 *99! ) 

991 FORMAT <* Y e 3.LE.yE2*) 

461 CONTINUE 

WRITE { 69 1004) XC?fDNUV*ONUC* TAH* TAUC* TftURV • TAURC 

1004 FORMAT* 32x#x c 2 =*£l2.4/32X*DNUV s *E 1 2 . 4/32X*DnUC 3*El2.4/ 

132X*TaU s#Ei? # 4/ 32X*TaUC =s*E12 .4/32X4TAIJRV »#Ei?*4/32X#TAURc = 
1*F12.a///) 

WRITE ( 6 9 1007) 

WRITER, 1010) HOL(l) >XV1*YV1 9 ZV1 , THV1 9 FMyl »DMUV tPvl t TVl * HOL ( 2 ) »XC1 v 
1 YC 1 * ?C 1 ?THC1 t EmC I ♦ONUCtPCl »TC1 *HOL (3) f X A * Y A , 2 A * THA 9 FM A 9 0NUa , PA 9 TA 
1 * MOL (4) *XB*YB»70* 

lTHB*EMB 9 r)NUC,PR,TB,H0L(5) . XC , YC v ZC » THC 9 EMC f ONUV , PC , TC , HOL { 6 ) vXD, yD 
1 9 ZD * ThU * £MD t 

1 rNUD * PD 9 TO 5 HOL { 7 ) * X V3 1 YV3« Z\/3 *thV3 »EMv3 9 DNUV ?PV3 ♦ TV3 *nOL { A) «Xc3»YC 

13*7C3,THC3* 

lFMC3 f nNUC*PC3 t TC3 • 

IF ( XF # GT^XFT) WRITE {69 1005) 

100^ FORMAT <2nX*NO F POINTS) 

TF 1 Xf.Le>XfT> WRITE 1010) hOL fo) 9 Xp t YF * f thF s £Mp « fllMi t n « p F , Tp 
• IF fXG.GT.XGT) WRITE (6,1006) 
lOOk FORMAT <2oX*lMO G POINTS) . 
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IF (XG.LE* XGT> WRITE (6,10IO)HOL (10) « XG , YG 9 ZG 9 THG t EMG t DnUO * PG * TG 
TF UTaU.fQ* 1 ) WRITf ( 6t 1010) HOL (1 1 > 9 XV4# YV4 9 ZV4# ?HV4t jtMv4 9 pNUV4fPV4 
1 * TV4 

WRITER, ]013) 

1013 fORMat<///) 

IF(ISh,EO.I) GO to 7384 
WRTTFf6f 1007) 

1007 FORMAT ( 1X*PT*<*6X#X# 11 X#Y* 11 X*Z*8 X#THETA#7x**mACH» 8 X*DNi |2D* 10 X # P tt UX 
1*T*/) 

WRITE (6f 1010) (HOLE (I > 9 X 2 (Ii ,Y2(I) 9 Z 2 <I) fTH2 a ) *EM?(l) *ONUE ( I ) 9 
IP* m .T2 (I > 1 1*1 .96) 

1010 FORMAT (lXA2f8El2»4) 

7384 CONTINUE 
Xl=Xvf 

CALL pOLY (XVI >YV1 tTHI *XV3«THV3, XV2 *THV2?XC1 • Ycl *XC?*T hC 2 ’ x 1 > 

CALL THRUST (XVI 9 YVl » EM VI * XV3 9 YV3 9 THV3 * X V2 ? X V4 9 T AU » o^Uv » XC 1 « YC 1 0 
1EMC1 , xC3,yC3 # XC2 9YCZ9TAUC.0NUC, iTAUtTHCl, THvi t THx> ' 

WRITE (6*09) 

WRITE (6f 99) 

GO TO 6000 
462 CONTINUE 

WRITE (6*992/ 

992 fORMAt(*THV2,LT»0*) 

450 CONTINUE 

6000 continue 

5000 CONTINUE 
END 
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SU 9 R 01 JTI £ CMA° 

c 0 MM OM / ft / P 1 * T m I 9 £ M i * G I ’ T I * H T. * •■'MH I * U l 
r OMMOM/^/ X H » Y A , X - • Y d * X c ♦ Yc * Xq * vp, ♦ 

1 Tp A 9 1 Hci * f MC * [rih t £Mjy ? 9 uMq « g-Mp. * X M U A f XMU^ 9 XmUq « XmLJq 

COMMON/ C/L)NUrff nf-4Cic*xCl 9 Yd »xVl * Y^ 1 iL 

— COMMON/Mr/P A * pn , pc 9 PO , T A , T8 9 TC , TO • 6 A 9 r,« 9 6 C 9 GD 9 UA *UB , UC 9 UD 9 .-’HA 9 
1 , R H 8 9 R H C 9 f< ii U 
r OMMOM/|,/ Phi 
■ -yMijl = , .SIM U */r_' ; X > 

THAsTHi ft KMAssRMI $ UAsUi 

OWlJA-, • 

ca LL FlX (XCl 9 Ycl «TH 1 9 X M UI 9 XV i 9 YV l 9 THI ,XMUI, 

1 y A 9 y A t f H A ? x MU a * £ A 9 G A 9 

1 p A * T A 9 P H A 9 U A 9 p k ; ( I A ) 

‘ T H q = T h A - r> M u (J 

Cr.LL pm (DMUC 9 Pp*Tb 9 «iMB 96 R 9 UB 9 PH 8 ) 

F m p 1 . 7; p M[J 

' ' y M ijH '] sA^IN l l , / F M R 1 ) 

CALL FUlXtUrcii jHiS 9 XmUuI 9 X A t V A 9 JH A 9 XMUA • Xn 9 YB9TMB9 x m ug»EmB*gs* 
1 PH, TA.KHR9 UH 9 U'-JUC) 

Thc-ThI + L) ,,J U V 

£ A L L p M ( p N U V 9 P p 9 T q 9 hi ‘ v i c 9 6 f 9 U f 9 P u, c ) 

F M CW M C* 

• yMLlC ! - a 0 T N ( 1 # /CMC 1 ) 

CALL F A X ( X « ? Y A , TP A * A MU A 9 A V J. 9 Y V l. * THC 9 XmUC 1 9 XC * YC 9 THC 9 XmUC » £mC 9 GC 9 

1 r: ’c ’ I c 9 h hc * ^ * jj KjU v ) 

6U ONOUan^Uv+UNUC 

CALL pH ( ijNUu 9 pn * TO 9 tM!) 9 GO « UO 9 PHD ) 

T HO- r H i + 0 NOV-“|)MUl 

c A i-L F.l X ( X fci 9 Y , TH h , X Mu i ;i 9 X c 9 Yc » THC 9 XMU C * X D 9 Yo * T H P * XMUD , EMO 9 Go * 

1 P[) 9 T : J* P| )Q « u t ; , [;- iUn ) 

R r* 7 U H M 

— Ci\iP 
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SUBROUTINE THRUST < XVI fYVltEMVUXV3tYV3fTHV3tXV2fXV4«T4U»DNUV* 
1XC1*YC1»FMC1 , xC3»YC3f xC2,YC2.TAUC f DNUC* I TAU, THC 1 , THV 1 f THX ) 

COMMOn/D.i/ A2*PZ?CZ 
COMMON /V/ AVt«V*c v 
COMMO^/P/ PINf 
cOMMOn/SHF/ xshf,yshf 
cOMMOm/POL/ IPR*IP0LY 
COMMOmZ TnRMAX/THPM ax 
cOMMOn/PmPM/ PmOm 
100 FORMAT 13*5) 

' GETZ tx5 ) rAZ*X5#X5*BZ*X5+C Z 
XL T Fen • 

XTH~0 » 
y M 0 M sr 0 « 

DX= (XvS-XVl) /5o* 

Ns(XV3^XVl>/0X*l. 

DXNs: (XV3»X VI) /FLOAT (N) 

N»M«- i 

CALL pM<DNUVfP M TA*EM A *<3 A tUA*RMA> 

XAcXvi $ Y A = Y V 1 % ZA=geTZ<XA) 

THftaTHVi 

CALL Fix (XA*YA* ®7854*0* *XA,YA9n* * 0 • * XA f YA *THA • XMtJA t EM A f GA ♦ PA f TA 

lRHA*Uft»DMUV) 

PAsPA»PTNF 

PX*PA*PJNF 

PX=PX/PTNF 

IF ( I PP o FO • 1 ) GO TO 30 
WRITF(6*1 10) 

11.0 FORMAT (///35X* VEHICLE UN£>ERSURfaCE*/7x*X* 12X*Y«12X*Z*9X*P/PINF* 
l A X*M Ac H*9X*LlFT*RX*THRUST*7X*MnMENT*) 

V/R I T E ( 6 * 10 0) xAtYA*ZA,PXtEMA*XLlF.XTH,XMOM 
50 CONTINUE 
DO 1 o 
XrsXa+DXM 

Yr=Yv1+ {XR-XV1 > / (XV3-XVI ) ^(YV3-YV1) 

7r=GFT/ <XB 5 

THr- A TAN( (YV 3-YVl)/(XV3-XVn ) 

CALL' PM ( DNUV ,pp f T8*EM0,GB*UB f RuB ) 

CALL FIX (XBfYBf •7B54»0» »XB*YR*o* * 0 • * Xr t Y8 « THB * XMUB * FMr #GR * PR • TB 
1 RHR * Up 9 r>HUV ) 

PBaRB-PTMF 

HUM* (pA + P8 J # (Zr + 7a) /A* 
nXLF= DUM# <Xb-XA) 

. cXTHsrUM^abS ( Y p-Y A ) 

Yrr = ( YA+YB ) /2 • « YSHF 
yBB=-{ (XA+XG) /P b -xSHF) 

OXMOMa-YBB*DXTH^XBB»DXLF 
XLIF=xLTF+QXLF • 

XTH=XTH+nXTH 

XMOMsxMOM+DXMOM 

X A - XB S Y A - YB * ?A = ZB $ PA^pp 

fma=emB 

PXsPa+PIMF 

PXaPX/P IMF 

IF(IPP*E0»1) GO TO 31 

WRITE f 6,1.0 0) XA*YA*ZA,PX,EMA f Xj.lF t XTH,XMOM 



Page A3- 11 


TR 188 

51 CONHmUE 
10 CONTINUE 

Ms (XV?“XV3) /DX*1 * 
nXN= (XV?-XV3) /fLOaT (N) 

NIbN+1 

no 20 I=?»N 

XR=Xa*DXM 

YRsAV+BV* (Xb«Xv3) + CV* (X0«XV3) 

70 kGEtZ { xB ) 

TH8 = AtAN (BV + 2**CV<* ( X B - X V 3 ) ) 

TAURRsTHR^TNV3 
RATB= (Xb~XV3) / (XV2-XV3) 

IF ( ITaU # EQ 9 1 > PATBs(XB-XV3)/(Xv4*XV3) 

IF (RATB.fiT*l* ) RaTB-N 
ONUB-DNUV + T aURr -►? p A T8 4> TAU 
Xl*XV3 $ x2sXV2 
TF(ITaU.FO,1) x2=XV4 

rata= ( <xi-xci > / (xa-xci ) -j, » ) / c (xi-xci > /Ua-xci ) *i ,) 
if (Bata *ge • l • > Rata = u 

TAIIXsATAN{RATA*TAN(TAU) ) 

d NUB*dNUvO"AURr*?p*TaUX 

CALL pM ( hNUb *Pr « Tb * EMQ *Gb »UR *Rm8 ) 

CALL Fix (XB^YB* *7854*0. * XB • Yp 9 0 , * 0 . ? Xr * YB * ThB * XMUB • fMr • 68 * PR » T0 * 
1RHB«Ur*DNUB) 

PB=PR*PINF 
OUM* (pA*P8 ) *<Zr*7a> 
rXLF“ OUM#(Xq-*x a ) 
r} XTH= r> Uivi *ABS(Yp-YA> 

YRR- (YA+Yr) / 2 * « Y SHF 
X RBk-( ( x A+X0>/?,-xSHF) 

OXMOMe-YBR*DXTH*XBB # DXLF 
XLIFsyLir+OXLF 
XTHwXrH+oXTH • 

XMOM=xMOm+DXMOM 
X a -* X B $ YA = V0 % 7a = 2b $ PA “PR 

pTMA^FM© 

PX=Pa+PTNF 

PXsPX/PIwF 

IF ( I PP ® EO a 1 ) GO TO 52 

WRITE (6,1 00) XAfYAtZA t Px f EMA t XLlP“tXTH,XMOM 

52 CONTjnUe 
20 CONTINUE 

XL V*Xl If 

XTMVry TW 
XMOMV/sXMOM 
WRTTF ( 6* 1 0 I ) 

1 0 1 FORM4t(/> 

XA = Xci $ YA = YCi s ZAS6ETZUA) 

XLIFsO. 

XTHso , 

XMQMsn# 

Ns (XC r~Xc1 > /DX*1 . 

OxN=<X C 3-Xpl ) /FLOAT (N) 

NsN+ 1 

CALL pM (0 NUC»Pa »Ta»EMa *Ga9Ua ?Rha ) 

THA=THCl 
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CALL fIX <XA*YA* *7854*0. *XA*YAi ft. *0* «XA*YA»TwA»XMUAtrMA»GA9PAtTA« 

IRHA, UAtClNUC) 

PA”PA-P T nf 
PX»P/uPJNF 
pX-PX/PlMF' 

IF ( IPB.EO. 1 > GO TO 53 
wR ITE(B*111) 

111 FOPMAT<///50X#rOWL* /7X*X* 12X*Y#1 £X*Z#9X*P/PlNF tt 

1BX*MAcH#$X#LIFT*RX*THRUST*7X*MomENT*> 

WRITE (6, 100) XA*YA*2A f PXfEMA t XLJF f XTH, XMOM 

53 cONTImUC 
DO 3o 
XB=Xa^DXM 

YpnYcl + (XB^XCl } / (XC3*XCl ) * < Y c3*YC1) 

7B=GET? (XB) 

THB e ATAN( ( Yq3-YC 1 .) / < XC3-XC 1 ) ) 

CALL pm (DNUc t Pq » Tb*EMb*Gb •Ub»Rhb > 

CALL Fix (X8*YB« • 7B5A»0, fXB*Yp«0. »0* tXn* YB»THBtXMUB*FMa*GB»PB*TBf 

IRHB, UB»DMUC) 

PB“PR«-R I MF 

DUM« (PA^PB^Zr^Za ) /*>• 
nXLF”«DUM<MXB»yA ) 

OXTH=nLM^ABS (Yp»YA> 

YBB- (yA*Y 0)/2 b «YSHF 
XBB=:.» j { x A * XB } / ?p«x5HF ) 

DXMOM--YRB*DXTh*XB8*DXLF 

XLIFsxLIF+OXLF 

XTH=iXtH^dXTH 

XMOMsxMOM^DXMOM 

Xa=XB $ YA=YB % 7a=ZB s PA=Pp 

FMAsEM^ 

PXePft+PlMF 

PXbPX/PINF 

TF(IPP*F0*1) GO TO 54 

WRITE ( 6 9 1 00) XAtYA*ZA t PXf EMA«XI_lF t xTH,XMOM 

54 CONTINUE 
30 CONTINUE 

N* (XC?-XC3>/DX + i . 

IF (*C?* EO.XC3) GO TO 60 
DXNr (X C ?«XL3) /FLOAT <N) 

N s N ♦ 1 

DO 4 0 I s ? 5 N 
y B = X A* Dx N 

•• YR=Yc^* (XB-XC3)/(XC2-XC3) # (YC2-YC3) 

7B»GFtZ (XB> 

PAT= (XB-XC3) / (XC?-XC3) 

DNUB = r»NUc*PAT*7 a #TaUc 

RATA=( (Xc3-XVi)/(XB-XVl)-l.)/r(XC3»XVi)/(Xc2-XVl)«l,) 

TAUYsaTAN (RATA^TAN (TAUC) ) 

nNUB»nNU0 2.#TAUY 

CALL pM (DNUB^PP, TB*EMB 9 GB»UB 9 RwB) 

THRsAtAN ( (YC2«YC3)/<XC2^XC3) ) 

CALL fIX <XB*Y0 * * 7B54, 0* *XB*YB?n, »0. « XR • YB t ThB « XM{JB «fMb « GB tPB »TB * 
IRHRtUptD NUB) 

PBsPR-PInF 

nUM= (pA^pe) * (ZR*?A> / a. 
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DXLFb J)UM* <X0-xA) 

D XTH s n UM * A fcj S ( Yra ** Y A ) 

YBB = (VA+VB) /2,-YSHF 
yB8=*( (XA4-XB) /?, ^ySHF) 

OXMOM»-YRB # DXTH*XB0*OXUF 

*lifsxlif*dxlf 

XTH=XrH+nXTH 

XMOM»yMOM*OXMOM 

XA = *B 5 YA^YB * ?A = Z B $ pa=pr 

fma=emB 

PX=Pa+PTNF 

px=px/p jnf 

IF ( I PP«FO. 1) GO TO 55 

WRITE <6, 100) XA*YA,2A f PX,EMA,X| IF, XTH,XMOM 
55 CONTTM^F 
40 CONTjMUf 
60 CONTINUE 

XTH = XjH4 XTHV 
XLIFkXLIF+XLV 
XMOMsxMOM+XMOMV 
WRITf'( 6*4?03) vL IF 

4203 FORMA T (///20x^TOTAL LIFT ~*E12 # 4) 

WRITE (6, 1008) xTH 

1008 format.(/2ox*thrust (from cowl and vehicle unde RSur face> =*e12«4> 

WRITE 4206) xMOM 

420 fo FORMAT (/?o x ^ m OmENT (FROM COWL aND VEHICLE UNDER5URF ACf ) =<> El2»4) 

IF (THx*LT*l*E+n9) GO TO 4201 
WRITE (6,4200) 

4200 FOPMAt(/?0x^ExTT PROFILE AND TwRyST AnD MOMENT (FROM sIDE w^^-LS) N 
inr Cai ClJLATED,*) 

GO T 0 **?. 1 2 

420 1 IF ( T h X • G T * 0 * ) GO TO 4205 
WRITp (6,420?) tH x 

420 4 FORMAT (/20x*THRUST (FROM SIDE wALLS) ~#E12 a 4) 

PMOM ) r 0 * 

WRITE (6,4211) PMOMl 

4211 FORMAT </?0X#MOMENT (FROM SIDE WALLS) = **£12,4) 

GO TO 4?i2 

4205 DTHOstHX-XTH 

WRITE (6,4204) pTHQ 
XTH=jhX 

PMOMl =PMOM-XMOM 
WRITE (6,421 1 ) pMOMl 

4212 CONTINUE 

PTHL= (THRMAX-XTH) /THRMAX^loo , 

WRITE ( 6 ? 1563) PTHL 

1563 FORMAT (/40**>P£RCENT THRUST LOSS s oF6-,3) 

RETURm 

END 
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SUBROUTINE: XTH(X2,Y2tTH2fDNUE*PlNFtTHS2t 
1 THX 1 ? XMASS1 *THX) 

C0MM0M/0B/AZ9BT*C7 
COMMQm/SHF/XSHF * ySHF 
COMMOn/PmPM/ PmOM 

c OMMOn/oNU3d/omU3d 

DIMENSION XNU(FO) f X (50) * Y ( 50 ) *7(50) *Em ( 50 } t P (50 ) * TH ( 5o ) t T (50) 
DIMENSION DNUEUo) * x 2(lO)»V2d n )»TH2( 1 0) »RHO(50> tU(50)*G(5o> 
GETZ<x5)«AZ*X5*X!WBZ*X5*CZ 
thx=o 9 

TF(AZ*9Z.EQ*0.) return 

INdEX*0 

DO in 1 = 1*5 
Jl»(I-l)#5*l 

TF(X.fQ,B) J2=?6 

no 20 Us JX *02 

PAT = F| 0AT(U-J1)/FL0AT(J2-J1) 

X NU( J) =DNUE (I) *RAT*(DNUE (Ul) ~DNUF ( I ) ) 
X(J>sv2<I)*RAT#<X2(I*l)«X2(D) 

Y<J) syH(t) + RAT*<Y2(Id>-Y2<I> ) 

TH (J) -TH? ( I ) ♦RAT# (TH2 (I ♦! ) -TH2 ( I) ) 

TAIL pM(XNUCJ) »P( j> »T( J) ,EM( J) (J) »U(J) tRHO(J) > 

CALL fIX(X(J)*Y(J) *. 7854, 0, ,x ( jj #Y ( J) ,0. tO, ,X < J> * Y (J) ,TH (J) »XMU* 
lFM(J)tG(j)*P(J)*T(J> 9 RHO ( J) 9 1J ( j ) tXNU ( j) ) 

XNU ( J) =DNU3d 
7 ( J) =f*ET7 ( X (U) ) 

ao continue 

10 CONTINUE 
Q6 CONTINUE 
T 3 1 

f 1 bRH 0 (I)*U(I)*SIN(THS 2 -TH(I) )/RlN<THS2)*Z(n 

rlB<R H 0( I> «U<I)*U(I> *SIN (THS2-TH<I> > *C 0S (Th(I) > * <P < I > -PlNp) *SIN (TH 
1«>)*7 <I)/SIN<tHS2> 

Wl n (RhO(I) «U(i)Hi(D *SIN(-THS?-th(I) > *Sl N ( TH ( I ) ) * <P ( I ) -PI NF) <*COS (TH 
1 52 ) )<>7<I)/SIN(THS2) 

XMASS?=0, 

TH X 2 - o * 

PMOMsO • 

DO 150 1=2^26 

F2-PHO (I >*U ( I) ttSlN ( THSa^TH { I) > /SIN (THS2) *Z ( j > 

f2s(P H 0 (I)*U(I)*U(I) #$IN(ThS?-THU ) > * rOS (T h ( I) ) ♦ < P ( I > «*P I NF ) *S IN < TH 
1S2) ) £ 7 ( I ) /SIN (THS 2 ) 

WJJ=(R M Q d ) *U ( I ) *U ( I ) *SIN {ThS?-TH ( I ) ) *SIN ( Th ( I ) ) * (P ( I ) -PINF ) ^COS (TH 
ls2) ) *7 ( I) /SIN (THS2) 

THx2 = tHx?9. b 5^ (FI *F2) *(Y ( I ) -y (T-D) 

XMASS?=XMAS52**5<ME1*E2) * (Y (H-Y (I-{ ) ) 

DUMls # 5o (Fl *( y n-1 ) -YSHF) *f 2* ( y < I ) -YShF) ) 

0UM2s # 5* { Wl# ( X ( I— 1 ) ~X5HF ) *W2 # ( V ( I ) -XShF ) ) 

PMOMspMQM* <DUMf«nUM2) * ( Y ( I) -Y ( T-l ) ) 

E 1 s F 2 5 F 1 -F2 % Wl=W2 
150 CONTImUF 

PATM=xMASS2/XMASSl 
rMS“R * TH- 1 

TF (ARs (EmS) .LT. l ,F»03)GO TO 46?5 
INDE X = I NDEX ♦ i 
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TF ( Index, GT. n go TO 4625 
no 1 7 n 1 = 1 *26 
FM1=FM(FM( I) *G( I) ) 

FM2=FmI/PATM 

ITM=] 

£MT=R4TM*EM {1) 

197 CONTimUf 

PMTsFW <FMT *G ( I ) ) 
pRM=(F M ?-FMT)/FM2 
IF ( ARS (FRM) .LI. 1 ,F*03) GO To l7f 
DUMDssi' • / RATM 

call pRROR (6» ITM,g;MTf eRM,dUHd?fMT1»ERm1) 

GO TO 197 
1 7 1 CONTINUE 

DUM = 5qRT ( (G(I)4i.)/(G(I)-i,) ) 
frTSOsqQRT (EMT*eMT-1 • ) 

FSQ=So R T(EM(I)«EM(I)«l e ) 

XNU2=oUM*ATAN (ftsq/dum) ~ATAN (ET$Q ) 

XNU1=oUm»ATAN (fSQ/OUM) -at AN (f50) 

XNU ( T ) “X NU < I ) ♦ x N U2«*XNUl 

CALL pM(xNU(I) ,P(I)»T(I)»EM(i),6(I)»U(I> »RHO(I) > 
170 c ONTim u e 
GO To «6 
4625 CONTINUE 

THX = ThX?/RATM-.THXi 

RETURN 

FNO 
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SUBROUTINE FIX (XI f Yl *TH1 , XMUl * V2* Y2 • Th 2 i XMU?, X3 * V3, THs , XMU3f EM3,G3 
l»P3*T3*PH3*U3?ONU3) 

C 0 M M 0 M / n 3 / a Z*rZ*cZ 
C0MM0N/0NU30/DNU3D 
TTM=o 
99 CONTINUE 

SL1«TaN<.tH1*XMij1 ) 

SL2 S T ftN { TH2*XMU? ) 

CALL ©EM (XI tYUSLltX2*Y2.SL2tX3*Y3 ) 

IF (X3.LT.O, ) WRITE ( 6 * 100) 

100 FORMAT^* NO X3 POINT IN FIX*) 

IF(X3,LT.0t) STOP 
RQ?0 =rH3^U3 
200 CONTINUE 

IF (X3.LT.0 • ) GO TO 99 

7 = ft Z#y3<*x3 + BZ*x3 + cZ 

P030=R02n/Z 
DNU3D=PNU3 
JTQ=1 
FRT=,001 
10 CONTINUE 

CALL pM(0NU3Df P3tT3fEM3*G3tU3»RH3) 

RQ3T»RH3*U3 

ERO= (P03T-R030) /R030 
IF ( ITq*GF« 1 0 ) FRTs.005 
Tf(aq^(ERO) •LT.FPT) GO TO 20 

CALL' F^ROR(20 f TT0,DNU3D,ER0, 1.^ , DNU3Di. *ERQl ) 

GO To In 
20 CONTINUE 

XMU3=aSim U./EM3) 
rL13s5L1/2*^«5*TaN (TH3+XMU3) 

IF (XmuI.FQ.O, ) SLI 3=SL1 
SL23 e sL2/2 # * # 5*TAN(TH3-XMU3) 

IF (XMU2.EO*0* i SL 23 = SL2 
CALL gEM()<1»y1* SL 13ix2*y2?SL?3,x3«y 3) 
iFflTM.EQtl) GO TO 300 
ITMsI 
GO TO 200 
300 CONTINUE 
RETURN 
END 
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SUBROUTINE THMfRHl »UI«PI 9 RINF»YV 1 tYCl*DYV»Xv 2 *THRMAXfDNUVeDNUC? XC 2 
1 ) 

COMMQN/DVAZf B 79 CZ 
Al=YVl-VCl 
Fl«RHt*UI*Al 
22 =A 7 #XV 2 #XV 2 +rZ*XV 2 *CZ 
ZSkAZ^-Xc^XCS + rZ^XCS^cZ 

A 2 = (Z 2 +ZS)/ 2 .oDYV 

TT = 1 

nNU»(nNUV^DNUC)*SORT( (Z2*Z5) /2. > 

10 c aLL pM(nNU»Pa,T?»EM 2 9 G 2 fU 2 fRH?> 

ER T =fF 2 -Fl)/Fl 

IF (ABS fERT)' .LT.l .E-O*) GO TO 2 * 

CALL frRR0R^99*lTfoNU*F^T t«9»f)N(ll9ERTl) 

GO To 1 0 

20 THIS (Pl-PINF*RHI*UI*UI) *A1 
Th?= [p 2 -piNf^Rh 2 *U 2 *U 2 > # a 2 

THRMAy = TH 2 -THl _ 


RETURN 

FNO 
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FUNCTION FM(XM,G) 

F X1 = *‘V U , ♦ „ ) /2.*XM**2) ( (rt+1 , ) /£*/ (6-1 . ) ) 

RFTUR- ' 

F n r< 



SUBROUTINE Ut>1 1 X A 9 Y A * 5L A , X B * Y« , SLB f X C , Y C ) 
XC= ( Yw-YA + SLA^x A-SLb*Xt$ ) / (SLA"SLB) 

Y C = Y A ^SI.. A <MXc-v A ) 

R (*; T U R N 
FT NO 
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SlJBRrHJTltoE ERKnR 
TTsini 

IF ( IT^UT. It?) GO T 0 12 
WRIT- < b * 1 3) 

13 pORi'ir T (#p:WWOR T^ST #) 

. . R T T (6«2<J‘> I 
2 u FORINT U5> 

STOP 

\2 rr(lT^T.2) 60 TO 1 * 

FRl=rR 
X 1 * X 

v ~ y ^ F 

J F ( X , F 0 * X 1 > X = X + * ! '2 
PfTURm 

’ -i^ >0 = X 1 -*.»! *(X-Xi ) / (EH-FR l ) 
p-P } = r P 
vlsx " 
vayD 
RF. TURv 
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SU8K0UI INt POLv {XVl*YVl>THI»XV? 9 THV3fK'/2*THv2t 
iXCj * YCi * <£*:» f tic.?.* XI ) 

C 0MM(JW/X/All f Alif A21 f A22 i A?3f a 31 tA^Zt A33 
COMMUW/POL/ 1 HP • T POL y 
n I M fe. MS 1 ON * < i Q ) * Y ( 1 o ) »TH(lo) 

IF ( IPnLY •F ,J » i > GO TO 20 

I T - 0 

X ( 1 ) = x ^ J 

X ( 2 ) - X Y I ♦ • 0 1 X v ? — X V 1 ) 

X(3)syV3 
X ( 4 ) -x 
TH ( 1 J - I HI 
TH<2)sl HV3 
TH (3 ) =ThV3 
TH(4 )a lHvZ 

TH< i ) =-TH ( i) S TH(2)»-TH(2) $ T H ( 3 > =“ T H < 3 ) $ THU)-~TH(4> 

y n > =ycj 

N i = 3 

wpnrf^,].o(!) 

100 pORMf-.T l ///3 oX«VfHIc l E C 0(>H f.) I n A t f;S*/ 

1 ] .5 x* X X * T 0*4 A # x <* 1 5 X *COOrO I NAT F S * ) 

3 (J On in I = 1*N 
y 2 = X < t ♦ l ) -X 1 
Yl aX ( I ) -X I 

As (TAW (T h U+l)] -TAN ( Th ( I ) ) ) /2„ / ( XZ-Xl ) 

RsTAn (lHll+1) ) -Z-*A*X2 
C = Y( l) - r#X1-a*X1*X1 
y(I + 1) «A*x^*X*i+«*x2 + C 
WRITE (6, 69) X (T > i X (1 + 1) tAttliC 
' 6 * pORMA T ( 1 p X D f 1U) 

10 CONTINUE 

TF (It, to. 1 I Go TO 2 0 
■r T = 1 

x a) a X cii A < 2 )» x n ) ♦.oi*{ x v2- x vn 

V(3)=X^2 
TH(2)aTHC2 
T H ( 3 ) -TH (2) 

TH (2) --Tm (2) .d TH(3)s-TH(3) 

Y ( 1 ) avV) 

Ms 2 

R 1 1 t f b 9 l 0 1 ) 

101 PORN- T ( /3()X*C(>nL COORDINATES*/ 

.1 1 ^ X * X #4 x # V U « 4 x X 4 1 S X *c 0 O W (j I M 4 T p $ * ) 

no TO J- 
2o rONTlMUE 
R F T U R \! 

FM0 
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SUBKOUnNfc PM (Ol-'li.UP** T2 t £M2 , G2 , m2* RM2 ) 

COMMON/ A / P 1 * T h I f E M i f G I » I' i « H I ♦ i: H i f U 1 

COMMON/ 0/ PHIi 

nTH=P./S7. J 

XFANsn^U/U I h 

T F < 0 T H • GE • U fN U ) !) T HsrL)N'J 

IF ( D i H f £ Q • u I’'* U ) TF A Nsi 

P 1 i=P i 

T i =T t - - 

FMIslmI 
HI =H.J 

" G ! =GT - 

pH j = rH1 

AisSopr ( gi*pi/phi > 

Ml =UT 
! ' 1 s fr" M 1 A 1 
1 1 1 = U 1 1 

mT=H1 +U) /d * 

P l “A L. At? ( P i ) 

XWU = r # 

• • X v l H g n * 

HO i 1 = 1 *1F Aim 

IF ( I IF Aim) l»TH = UNU-XNU 

XMU] =,,SlM ( i e /j£M1 ) 

R l = 0 ) *tMl/COS ( y NM1 ) 

P Z = - R ] w f > T M * P 1 
RH2 = ( P2-P i ) /G I 
RHZsrtH i *F.xP (PHp ) 

P ? P = “ X P ( P Ff ) 

PJ P=.7yP ( PI ) 

i |2 = Hl „d (t?j -I a ) * (P 2 P/HH 2 ^PiP/RHl ) 

! I 2 r: H T « u R / tl a 

T? = f- T (P?P»PMI 1 ,H?) 

G?“F oAH ( r 'd *H<dP* PHI 1 ) 

A 2 = (J?»P 2 P/HH 2 
F M ? = Sri*T ( U c/ a 2 ) 
y NUaXrviUl +L) I h 
Xm"1,=v K( I) 

Pl*P? 

Tl=T2 

01=0 

111 = U F 
PHI =rh2 

10 com 1 TmUi : 

dR = K >. p ( p ' ) 

M2 = S pi { I j 2 ) 

R IT T U R h ' 

FrO 
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fu^ct t.uu fi ( P i , f • H5 ) 

DATA t 63/*V 
T F L /v 0 s 9 

P = P1<M .0l3^SE + 0i>/2H6. 

H=H5 /iu . ff>3V/l.E + Q6 

F2»F*f 

A 1 ()=A|_.Un (R ) /2 # 3-5. 

79 = . ] ? b #■ A 1 U«Ai o -.27B*AlO 

?T=1 ' 

TF ( Ibl.liO, 1 ) GO TO 10 0 0 
T 63=: 

"■ T S 1 5 P n * ' 

T 0 = 1 S 0 0 . 

IF(F,fit.o#) Gu ro 120 

r = 60;i „ 

TOtrT 

X 0 0 u r QMT T kjUE 

IF (F.lT. i . ) <30 TO 4oO 
^0 To 1?. 

5 0 F 0 = ( H_til ) /rl 

I F ( A B S ( E ) • L T • l . F ~ 0 4 } GO TO 34- 

50 0 T sTn*l.i 

502 T T = 2 

IF (F,lT. a* ) OQ TO 400 
r;0 To l?n 
loo ns(^jin/^ 

IF ( A 3 $ ( E 1 ) * L T . I • E - 0 4 ) GO TO 34 
ITslM 

IF ( IT ^IT /d l ) oO TQ 10 
- - if (AR$ (il^uoo, ) ,LT. L D. ) GO TO n 30 
wR I T - ( o , f, 3 i ) Pif hS ♦ T 

- y 3 A FORMAT ( * LKkOK ID F I */# Hi = *r- 1 3 , 5 * 5v * *H l _ *E 13 „ 5* b\ * # T = *El3.5 
1 > 

STOP 

03u IF < IFpAG.tw. 1 } Go TO 5 U 4 
I f» L- / 4 R s i 
TO -200 0 ^ 

T — 2 0 0 0 • 

lFl^ M 'I.ij.) Go To *00 

r; 0 T !.,* A ? ■ 

5-o 4 TT ; '(0»li^ 1 1 

11 f -OPh..-.TU 1 t riP t Ri ‘ TUK t I N I T s £ ] T0 2000 - E RROp = * F 1 3 » 5 T 

GO To 34.- 

10 T^a T 0 -t> 0 <K 1 - r„ ) / ( E 1 -E 0 > 

jrOsj/ 1 

T 0 = T 

t=tr 

TFfF.L T* r.) Gu TO 400 “ 

120 Asl ,f_U7» (-. ii)42*F2 + .H242*F+^987> 

R = , 0 0 1 ** ( s 0 i i G / <* F ? ♦ * 1 5 0 3 *f H + . 9 3 r ) 

£ s - « r ?.y a v f 2 + • g / 3 1 ■» (- ♦ , 4 2 ^ 3 
IF ( F . I t , J. , ) Go TO 190 
As i . E- „7# ( 1 • j -5*4 
R=.0-)] * (-. ipG/^|: 2*1.1 1 [ t ) 

C=« „2 a n-iH 

190 TF U ,| t. 2000, ) ‘iiO TO 290 



As.OuoOOt* U.7 92#P2 + ,39b3«F+ } ) 

0= • 0 0 ] . * ( -9 • 0 5#p2 - . 07 9 1 7*f + . 245 ) 

-• Cs;lO,RO#f2 1 1 83«P+*97 

IF C F . » t . 1 . > Go in 2*0 

A=*000UU1<M4. 6 i -13,9»F ♦! 1 .c.q) 

-••- ns.CKil* (-2 j.ObtfF? + b6,tf 2iH"-52 *M ) 
C = 27,rt> # F*' "73* /3<>F*58.39 
C9° HlBA»T*T*b*l*C 
" ‘ T F n * L t • 2 0 U 0 • ) 0 0 TO 370 

Hlah1#a.+ (U+F)#(T/2000 p «l.)«?9) 
3 7 u CONTImUF 

00 TO ■■■ ■ — 

4()U 7 2=T*T 
T 3 = T 2 ^ I 
‘ T 4= T 3-s* r 

T5*T4 « T 

IF ( F , | i « - 1 * b ) no TO 45 0 

' ' !< M M 1 s 1 # 0 4 *i 

A1 s 4 , ? 4 9 7 6 7 B 
A2 = -& .yi 26 tJb 2 i £-03 
A3 s 3plbF??lJ4e-05 
A 4 = “ ^ . y 7 1 5432 £- Oh 
A 5=9.FU>3bttOE-i 2 
A6=-l r 0 l^6b32t+ 04 
F 0 T 0 4 1 1 
*50 CONTIMUF 

Alai „ 1 2 r,2456 
A ■?. e 1 , 3 9 0 5 7 1 0 E - 0 2 
A 3 s 2 s F b 6 b 3 7 4 ii - 1 6 
A 4s®] , 1 5 6 U 2 7 2 F «* 0 8 
A5sb.? J B fe9 c 9E-1^ 
a 6 a 5 ^ 3J;:P8Vbt + / ; 3 
Vmm1s?«, .54 

iybs) M 1 =A T ! + &2*1 * ♦ A 3*1 3/3. ♦ A44T4/4 

H 1 Sh .1 iiH 3 1 4 • / A-* 1 M 1 / 1 * t + 0 0 
350 IF (IT to 1) GU TO 50 
'~0 T(i* 100 
34o T (i ” T 

"■ fT = T#1 * r 

Mr T U * !- 
p'Ajn 
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F U M C T T ow RhEQ [ H , 9 1 * F , J ) 

T 1 =F T <F1 . F»h) 

T = T 1 / 1 • R 

P = P1 tt} ' Ql32S£+ 16. 

X F ( P f \ *' a < • e ) 0 U 'in 226 0 

KWM= i.b3*F*F-t> +78 .965 " 

fr N (Nt = j -if.) # 6 <5- f + 3 3 ,6 

IF ( T.pT ,2000 . ) GO TO 2Q3o 

" X H = F >\i M 

I F ( H . j I » 1 ► ) G vj To £16 0 
XM-FON 

■ GO TO *.160 

2 0 3 0 FF=F»f 

A = - 2 , 3*F F + 4 . i) 1 + 1 . 7 36 

~ 8-8.6] *FT-lb,4j?#F..6,66 

0 = - 1 6 . f* + Gs5 . 2 J . #F + 14.66 

X|M = - * 4 J7F» tt F F + * ; '625 tt F + 2 . 06 
D = A # ( A L 0 G ( 6 ) / 2 « 3 ) 1 <* 5 + B ( A L 0 6 < p } / 2 * 3 ) 4 C 

X.MrF fvH-0 «■ ( ( T-2 o 0 0*)/l00D»)**XN 
T F < F . L T • ) • J OO TO 2160 

As-,,«?i; w FF + 2.3F3 fr F* 1 .905 

fi* = ?#7ft*FF- / .bt>*F -6 .68 
ra-J. ft *FF+ / *3 0 «-F + 2/ * Uj 
X M = * « 4 7 # F F + 1 * 6 0 6 * p + * 3 5 

!) S A« (ALGO (K* /2,3 ) **i .5 + 6* (AUOG f P) /*,3) *C 
XM=FWW-h* ( (f-*000. ) /100D. ) *#XN 

r.O TO 2160 " ’ 

2260 KF=F^b 

1 F ( Kp ‘ fe 0 . - 1 ) \Ms 16.043 
TF (KF.tP.-2) XH=2 h. 0b4 
216 0 PHFQeP * x (•’ / T/d J 1 4.3*6. 2428 1 -02/32.1 74 
T a: T * * # 6 

P f TUR' 

r- no 



FUNCTION FtrAM (Tl »Pl *F) 

T=TX/} .8 

T ? a 1 # I 

P = P1*3 .03 32h£ 4 r,S/?l i6 # 

Vm = 0 . 

I F ( P . | T # 0 • ) tiu TO 5 d0 

T F ( T • !_t • 1 1)00 * ) ^ TO 4-4 q 

X ^ „ 1 5 p; - U B ^ I ? + • 0 0 0 0 9 1 tt T - » 0 b u 5 
44 0 xWs^F-n^J 2 OOOOei^T-. u 19 

IF ( F , !.. t , 1 • ) Go TO 4/0 
X Fjs * a 3 3 g $ U R T | T ) *“ * 0 0 0 3 9 1 * T * * £ 8 1 
4 /u G = «,l .A 33 £ *0 7*1 ? 4 # 0 U 0 o 7 £>#T 4 1 , 3 (R 7 

1 F ( T . 1 • ROO • } TiU TO 5?o 

G b 2 o £ ~ 0 j. s * T 2 - . 0 0 0 I 3 tf * T + 1 1 4 2 3 

.T F { T j * .2000.) GOTO 5 2 0 

Gb7. 5g /£-o^ # T2 -• 00 u4S7*t + 1 • B5 
520 0=0+ XM* ( aLOg (P) /?, 3-5* ) + XN* (p-i „ ) 

GO Tn b 3 o ' -■ • 

bbu t3sT2^T 

T 4 s T ^ 4 7 

C.P-Al + A£wT + A 3#T2 4 A4<M'3 4 A5#T4 
G»CP# (Cp.i , ) 

530 cONTjnU^ 

FGAM* G •' ' 

RF TOR? t 
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FUMa JOM Rl(Pi,K*Tl) 

Pa PI <'<■} * o 3 32bj.. ♦nb/2 116. 

■ ‘TsTi/i.n 
F£=F *F 

if (k,i.J . ) go Tf) ^oo 

t F ( T n I ,?0i)Q # ) GO TO 190 

T F ( F ; r, 1 . , , ) byj 7 0 

i2u a b I»E-'J7« 10-'*2*f 2 +,82 / ^*F + ..q87) 

Pa.(U«l*( # yU67»F2 + .1S«)3*F-> # 93 k> 

C = -* «F 2 +*6731*F + , 4 P.93 

GO T 0 £ Q ;' 

1 9 i Asl.F-C'7* (1 ,7G7<*r 2 -5 „46*F+5 .4) 
fj s # 0 (j J * ( - * 1 H 6 f * F ? +1 s I 1 # ^ + . 1 7 6 } 

£=-♦ + 3 # 9?5<*F -2*808 
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